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| EEE gemäßigt turbulente Bewegung. 
A Kae (Abhandlungen zur Hydrodynamik XII.) 
ER Von @eorg Hamel, z. Z. in Landshut (Bayern). 


Bei einer ebenen, stationären, gemäßigt turbulenten Bewegung in einem Parallelstreifen gelingt durch 
Einführung der Stromfunktion n) "neben « als Unabhängiger unter Vernachlässigung der Glieder mit v? die 
hypothesenfreie Einführung einer Turbulenzfunktion T?(y) — y das über x bei konstantem y gemittelte y —, 

von der man weiß, daß sie am Rande 1 ist und sonst >1, es sei denn die Strömung laminar. Turbulenz- 
funktion und. Verteilung der Geschwindigkeit berechnen sich gegenseitig eindeutig. Das über x bei kon- 
stantem y gemittelte Druckgefälle erhält keinen Teil von der Beschleunigung her. 


ER Regarding a plane stationary moderately turbulent flow in a parallel strip one succeeds in introducing 
Ba a turbulence function T? (Y) without any hypothesis by the introduction of a function y in addition to x as 
independent variable and the negligence of the terms in v?. Itis known that T? (5) — Y meaning the average 

of y computed with constant value of y — has the constant boundary value of 1 and elsewhere values that 


ERROR 


valently computable from one another. No fraction of the mean value of the decline of pressure, - 
with x varying | and y being constant, originates from the acceleration. 


DB: cIYyAe MIIOCKOTO, YCTAHOBHBITETOCH, "ymepeHHo BUXPEeBOTO BIKEHTSTMERA HaAPANeIb- 


HBIMH INIOCKOCTAMU YAAeTcH, BBONA HApANy C HesaBucuMoä mepeMmeHHOH x PYHKULIO TORA y, 


Ei; . -  H mpeHeÖperas WIeHAMH C v2, BBECTUH HECBA3SAHHYPW C AOHONHHTENIBHBIMH THIIOTE3AMH PYHRIHIO 
—- Buaxpa 7?(7), re y — cpenHee sHayenme IepeMeHHoÄ y, B3ATOe IIO x IIPH IIOCTOAHHOM y. 
OTHOCHTEeNBHO PyHKIHNH T WSBeCTHO, YTO OHA Ha TpaHnue pasHa ], a B apyrux Toyrax >18 
"cAyyae BUXPEBOTO IIOTOKRA. Dyurıns BHXpPA u päacnpepeleHue CKOPOCTH BEIYHCHAIOTCH 
B3AHMHO ONHO3HAYHO. CÜpenHee 3HayeHHe TaNeHHA NABJIeHNA, BEATOE IIO x IIPH IIOCTOAHHOM Y, 

He 3ABUCHT OT YCKkOpeHun. 


I. Wir leiten eine strenge mathematische Folgerung aus den Stokesschen Gleichungen 
ab Ar Hinzunahme neuer physikalischer Hypothesen. Unsere Voraussetzungen sind: 


1. Ebene stationäre Bewegung einer inkompressiblen zähen Flüssigkeit im Gebiete 
— o<xz<+mw, — h<zy<sh (Gleichungen (1) oder (3)). 
2. Die Geschwindigkeiten v,, v, bleiben beschränkt, ebenso Ay auf jeder Stromlinie. 


er für jede Stromlinie y = const. beschränkt. Das 


nennen wir: „gemäßigt turbulent“. Wir können deshalb y und x als Unabhängige einführen, 
y=F(«x,y) als Abhängige. 
4. Setzen wir Ay=fly) rw (y), 


so hebt sich aus der Bewegungsgleichung für die Stromfunktion y das », die kinematische Zähig- 
keit, einmal heraus. Mit dem dann noch bleibenden » gehen wir zur Grenze v—0 über. Anders 
ausgedrückt: wir behalten Glieder erster- Ordnung in » bei, vernachlässigen Glieder höherer 
Ordnung. Wir erhalten so aus (3) die Gleichung I, zu der II als umgerechnete Definitions- 
gleichung (4) hinzutritt. Rechnen wir (1) um, so entstehen die Gleichungen III; eliminieren 
wir aus ihnen den Druck, so entsteht wieder I. Die Ableitungen nach x sind bei konstantem y 
gemeint. 
Wir arbeiten nun mit der Mittelbildung 


3. v, ist im Inneren nie Null und 


[4 


ee | RL O0” 
1520 
—l 


wobei die Integration auch bei konstantem y gemeint sei. Insbesondere spielt weiter das so 
gemittelte y eine Rolle 


a) und TERN 


Su 


.. . are >1, unless the flow is laminar. The turbulence function and the distribution of the velocity are uni-. 
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In vollständig strenger Weise folgen jetzt die Sätze: 
AM (= 2) ist von y unabhängig, also konstant. Es sei positiv. 
0 . 


0% 
SER FE Boy re es: "(c), 
wo By=M & ne ist. 
Yy 


Bemerkenswert ist, daß das bei konstantem y nach x genommene mittlere Druckgefälle bei 
»—0 Null wird und keinen Bestandteil aus der Beschleunigung erhält. Das würde auch gelten, 
wenn die Vernachlässigung höherer Potenzen von » nicht gemacht würde. 


de BHIE AT. ee Dr A (d) 
Elimination von f und Einführung von 
dy=A(y)dy ; 
gibt die Gleichung 
=. v = —AB er ee re EEE (e’). 
v dy? 
Die Schwarzsche Ungleichheit gestattet nun den Schluß 
"det 
Das Gleichheitszeichen gilt nur für die laminare Bewegung. Wir setzen 
AB=T? (9) 


und sehen diese Funktion als Maß der Turbulenz an. Man kann noch beweisen, daß sie an den 
Grenzen y=y= Äh gleich 1 ist. Sonst ist sie unbekannt. 
Aus ihr bestimmt sich die Geschwindigkeitsverteilung nach IV, wenn wir noch Symmetrie 
voraussetzen. Umgekehrt bestimmt sich 7? aus gegebenem A durch die Differentialgleichung V 
mit der Randbedingung 7?=1für y= + k. 


II. In dem ebenen Parallelstreifen 
—o<rı<+mw,—hsys<sh 


mögen wir eine stationäre Strömung einer inkompressiblen zähen Flüssigkeit haben, die den 
Gleichungen 


9%, ov, ep 
A tee dv, Sn 

apa EA Zee 1) 
09V, 99, °» ( 
PIE Te en 


gehorcht. Es seien v,, v, die Geschwindigkeitskomponenten, p der Druck, o die Dichte, » der 
Koeffizient der kinematischen Zähigkeit, A der Laplacesche Operator. Führen wir der 
Inkompressibilität entsprechend die Stromfunktion durch 


9y Oy 
sry’ Uni PS De A ee 412) 


ein, so ergibt Elimination des Druckes aus (1) bekanntlich 


oAy Y oAy day 


ey ITS tes AA. oe are ee (3). 
Bei v»=0 wäre Ay=f(y). Das veranlaßt uns, 
Ay=f(y) +rw(x,y) 
zu setzen, indem wir x und y als Unabhängige, 
Yarılay) or m Re (5) 
als Abhängige einführen. Das setzt voraus, daß im Inneren = e nirgend Null ist. Eine 


are 0 
solche Strömung, bei der außerdem noch längs jeder Stromlinie d y/dx beschränkt bleibt, wollen 
wir „gemäßigt turbulent‘ nennen. 


Wir setzen weiter voraus, daß » klein, d.h. die Reynoldssche Zahl groß sei. Demnach 


_FauFr —2H,,F, F, aut 
Be EZ 


| i=- re FR, Er a +2 
i x 2 Fy 


"nursannnen.uen =) Hr Y- 


£ d En 
dw dw aa oy dw Oyoy_ | ” 4 I ] 
Te, är dl dy ap ay er ui) 
i FREE a 
eo. EEE RER 
r A, ie F, 
 vist herausgefallen. Setzen wir 
; } 1 en , l 
Ri 
so ist 
00_2F,F,, 1+F, 
Rn et ee, 
_ Deshalb läßt sich (4°) auch : so BR, 
Er | Rust Re ET Me 


III. Ehe wir weitere Überlegungen ve die Gleichungen I, II anstellen, wollen wir auch ge 
Gleichungen (1) unrechnen. Wegen 


a 
9%) y=const \0%) y=const oy 9%) y= const oy) F, 


ee, 
OYy) z= const oy oy F 


elle er ; : : 
_ werden aus (1, wenn wir hinfort unter 5, /jmmer die Ableitung bei 


EEE BEN 5* 


’ * 


air == 
er Fi 


. 
- D 


En ee (da 


iininert man hieraus ?, so a man 


A ER ES Re is Er] 
Ind FE ar wrf 
oder 
D:TH,Fn a un ‚1+fi 5 aF,F,v 1+F% u) 
| r Be 5 Ba ar =»f' F, +rf F, wis yYoz 
oder nach (8) 
aA „L+F F 
Fr) vs = —p +rfF,Ay. 


Das ist natürlich bei Vernachlässigung von »2 2 I identisch. 
Man kann die erste Gleichung III auch als Verallgemeinerung der Ber- 
noullischen Gleichung 


REN Ur ) Zn 
2x F +3 (va +) == vf @ 
schreiben. Die Ableitung links ist bei konstantem y zu bilden. 


IV. Nun arbeiten wir weiter mit der Mittelbildun g 


= lim — JS - 
>26 
wo die Integration nach x bei konstantem Y Vera sei (statt bei konstantem y wiein meiner 
Abhandlung über die Streifenmethode in den Abh. der Preuß. Akad. d. Wiss. 1943). 
Naturgemäß sehen wir v,, v, als beschränkt an; yist es sowieso. Dann folgt sofort 


In Mo 0 Br MF= 3 


Br UFERTTEE 


Dw., 
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a re VE EEE SER EEE Ta Ba 


Also folgt aus der zweiten Gleichung III 
| op 2 | 
M— - = oder BE ES ER ae 
| yo a) | u 
d.h.der gemittelte Druck ist von p unabhängig. 
Wir können also 


als konstante s mittleres Druckgefälle einführen. Hinge es von y ab, so müßte _ p/o mit 
Y 


x unendlich werden, was dem Vorhergehenden widerspricht. 
Aus der ersten Gleichung III erhalten wir dann mit 


1+F3 
M IM EBD N ee Eee 9 
4) -2@=2W @) 
Er Br (e) 
da 
a 
Er er ist. 


Es istsehr bemerkenswert, daß bei der Einführung der Un- 
abhängigenyundxkeine Einwirkung der Trägheitsglieder auf 
das mittlere Druckgefälle $ stattfindet. Das gilt auch, wenn 
wir die Glieder mit »® nicht vernachlässigen. 


0 
Aus I finden wir nun, da M 5 = () angenommen werden darf, falls nur w auf jeder 


Stromlinie beschränkt bleibt, 
a 2 EEE ee ee (d) 


wenn wir 
Be ee ee ee (10) 
setzen. Gleichung II liefert, da M (G) = Bist, 
dB | 
—-— ME re Er 
JAH HM wR,) (@) 
d 
wenn wir M (F,.)= 0, d.h. F, = Er längs einer jeden Stromlinie beschränkt voraussetzen. 


Nun folgt aus (c) 


ge, 2 dBläy 
f = vB ’ r E7% y Ba ’ 
somit aus (d) x 
L:=B B 
A m — —— ı —, 
f a 7 
Verglichen mit (e), verlangt dies 
vM (wF,) =0. 


Es muß dahingestellt bleiben, ob dies die Vernachlässigung des Gliedes mit » in II erfordert 


oder M (wF,) =. 
Wir führen nun zweckmäßig ein gemitteltes y als Unabhängige ein, d.h. 


Per wP 
Ta = lim — a re aM TR le er ee Na 11 
dy= Ady lim | My de dy lim — j dydx (11) 
—i 4 


eben das mittlere y eines zu y gehörenden Teilchens. 
Damit wird aus dem veränderten (e), d.h. 


dB 
re: ER ee Mae Wen Viral Ana ie Bela a, et Ar (e) 
dB ’ 
en d7 et Re ee a re (e‘) 
ee S Mal: une _ ie (c’) 
gr Fr ee ©; 
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Hamel, 


Eine gemäßigt kubalonie Bewegung 


Über die Bedeutung von A läßt en folgendes sagen: Es en = 


2 d.h. die Ableitung der Stromfunktion nach ER 
Y . 5 & 


gemittelteny, ferner ist 


erstens — 
Fr 1 
A=M(F)=M 2”) - u(), 


0) Y z=const 


AUSH 188 a auch das Reziproke der gemittelten reziproken Ge- 


schwindigkeit v,. 


V. Wir haben nur eine Gleichung zwischen A und B. Aber wir können durch Benutzung 
der Schwarzschen Ungleichheit noch eine wichtige Folgerung ziehen. Es ist 
” 2 


({ nr) -| 2 Ba En ee 


nach Schwarz 


= er [re 
ii y 


Also 
ae. 
AB> (lim; | yi +F% 22) = (tim; = >; 
—u 
wo s die Bogenlänge der Stromlinie von = —.a bis x +a ist. Wir setzen deshalb 


AB=-TY)>1 


und sehen T?als Maß der Turbulenz an, nennen es „Turbulenzfunktion“. 
Für die Laminarströmung und nur für diese ist sie überall gleich 1, sonst größer als 1. Denn _ 


für die Laminarströmung ist erstens s — 2a, zweitens F, = = von x unabhängig und F,=(, 
weshalb dann das Gleichheitszeichen gilt. 
Auch am Rande ist ?=1. 
Denn zweifellos ist am Rande, d.h. für y=-+h,y=0 oder y=@ (Durchflußmenge), 
v, = 0 also r —=0,also B=0. A dagegen ist dort unendlich. 
Sei nun an einem Rande 


yzala)y+h’+Pp)(y+h’+. 


so ist Ay=2a(e) +H3ß(y+h)ta’(y+h) Heß +hj’ +. 
Ady=2u" +3" (y A) 2RT +5" (y =, 
Ay ay AAy Ay i 
Da aber für y= — h die linke Seite — — —— — Null wird, muß 40” = 0, a’ linear 
dx °%y oy 0x 


sein. Es muß aber sogar konstant sein, denn andernfalls bliebe v, = = nicht beschränkt. 


Ist aber in der Nähe von y= —h 
y=oa(y-h)®-+.... mit konstantem o, 
so ist 

a 
0x — 

Oy — 1 — 1 1 

— zw2oa(y-+h 97 re 4A:= lim — ETW rn [pa 

5 (y+h+ Yay+ De [: lat: 


ae | ° _— 
B=lim Sht 2222-21 +. 
das Produkt AB=1-+...,d.h. in der Grenze 1w. z.b.w. 


 AB= a * ergibt sich schließlich die Ba, indigkeitsverteilung 


en ehr) 


SR, "einer Differentialgleichung, die mit m — R für 3 %= + hintegriert werden EE Gese En win- 
 digkeitsverteilung und Turbulenzfunktion bestimmen sich ge- Re 
"genseitig eindeutig: RE 

Bemerkenswert ist noch folgende Formel: für die gesamte durchströ- h 
Fr mende Benes RR Aa 


d: B 


ame. 
EN RR NER 
an 2 
dy 2 
a a re VE a 
ist | * 
VI. Ist beispielsweise A= const, so folgt - 
Er: ER SA“ 
= | aan 
daraus a SEE 5 
MER "%, Br (© 2) = 34, 72 Pi 
7: RN ne | 
r T® ist das Maximum von 7°, bei Symmetrie bei y= 0 gelegen. n 
A -s Weiter folgt 7 Er 4, Zone s 3 a 
2 | = | 
3 = dr: 2 
£ == 
. 2 cv p®8 
i ; Frm: 


"Um ein weiteres Beispiel "zu geben, versuchen wir den An 
Parametern aundn _ 


1er jo<acı, a | 


Damit ist N | en 2 
- dy Y Fe, Yy 2n +1 2 a * £ jr v 
[#- ala el ) ; Ben 
8 (1a) (y Y ante, Aus 
ee 7 
‘0 - j we, 
ea filnae (1-0) = & SE 
Ber, | =2.| 3 Fear) On) Tonrıp (in 34 ee 
3 woraus sich die mittlere Geschwindigkeit zu nn 
= 
na Er °, 1—o)? a(l—a). a2 ) 
e Mn I iurss & om) (2n-+3) Fonrıy (dn+3 
x berechnet. Die Ge mekelteveritss berechnet sich zu 
3 dy_ 1 3m n ee > e Rn E 
3 FR ETEERN 1a Ei 4 ) 2 ( N) Temrn) | -4 ) B u 3 
Daher ist die Geschwindigkeit in der Mitte y= u | S | 2 
N 1—a & | 2 
| TR Mia), Tan) An >31 
| und das Verhältnis _ 
1—a & 
IE da) TEST 
Um ee! (10) a2 
+2. FI) ün+3) T On Fiwinrs) 


Sieht man dieses Verhältnis aus & Erfahrung als gegeben an, so berechnet sich, wie man leicht 
nachweist, für jedes n und jedes > — innerhalb gewisser Grenzen eindeutig ein & zwischen 0 


und 1, das mit n— oo gegen 1 ns 
Das gegebene Beispiel soll aber nur zur Illustration dienen; das wahre 7%? wäre aus A 
in der angegebenen Weise zu berechnen. Das muß ich aber mangels Unterlagen anderen über- 


lassen. In der Praxis wird man wohl y mit y und daher z mit dem gemittelten v, identifizieren 


dürfen. Doch bleibt das experimentell zu untersuchen. 
Eingegangen den 22.8.47. 
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. gen in einer unendlich kleinen Umgebung des Verzwei- 


- Zur Theorie ins an Verdichtungstoßen 
Von Walter Wuest in Göttingen. 


Es wird ein vereinfachtes zeichmerisches Verfahren zur Berechnung von 1 gegabelten Verdichlungsstößen. 
erläutert. Bei Annäherung an gewisse Grenzlagen, in denen einer der Stöße in eine Machsche Welle 
entartet, versagt dieses Verfahren. Die Grenzfälle werden daher analytisch behandelt, und die Gleichungen 
der Grenzlagen werden angegeben. Der Sri von Gabellösungen wird an ie verschiedener 

: Diagramme erläutert. 


tert A simplified graphic method of computing bifurcate Dee of compression is demonstrated. Approaching 
certain boundary positions meaning that one of the shocks degenerates into a Mach wave, the method is 


Jailing. The boundary cases are therefore treated analytically and the equations of the boundary positions 
er stated. The region within lach bifurcate solutions are existing is demonstrated by aid f different 
jagrams, % 


B arof pabote onnchBaeTtca rpadhnaeckuit MeTon pacyeta PasBeTBNeHHEIK Tunpapımyeckux 
yAapOB. ITOT MeTON HeNPHTONeH IPH IIPHÖNHSKEHHH K UCBECTHEIM TPAHHUYHBIM TIONOKCHHAM, 


KoTNa OnuH #3 ynapoB nepexonut B Boumy Maxa. Ilosromy 9Tu rpaunyHse cıyyan 
HCCHENYIOTCH AHANATHYHUCKH, IIPHYEM 1AI0TCH YPaBHeHHA IA TPAHHYHLIX NONOMeHnK. O6nacTk 
HPHMeHeHHS- peımeHusa 3aay PasBeTBJIeHHBIX YNAPOB HOACHAETCH ua HOMOLLH Paauanez 
Auarpamm. 

1- Einführung. 


- Auf Schlierenaufnahmen hat man beobachtet, daß Verdichtungsstöße, die mit einer Ab- 
lösung der Grenzschicht verbunden sind, sich im Inneren des Strömungsfeldes verzweigen, wobei 
der eine Zweig wieder zur Wand zurückläuft. Die Theorie dieses gegabelten Verdichtungsstoßes 


“ist von A. Weise[l]und H.Eg gi nk [2] behandelt worden. Auf Veranlassung von Prof. 


Tollmien wurde vom Verfasser eine neue Darstellung der Theorie des gegabelten Verdich- 
tungsstoßes gegeben, die gegenüber den früheren einfacher und leichter auswertbar ist und die 
Grenzlagen des Stoßes vollständiger behandelt [3]. 

Eine vollständige Theorie des Verzweigungsstoßes müßte die Eigenschaften der Grenz- 
schicht mit enthalten, weil offenbar Grenzschichtvorgänge an der Auslösung des Stoßes wesent- 
lich beteiligt sind. In dieser allgemeinen Form ist die 
Aufgabe jedoch nicht Gegenstand der Theorie gewesen, i MG PR 
die sich vielmehr darauf beschränkt, die Stoßbedingun- 


gungspunktes zu berechnen. Die im allgemeinen Ver- 
lauf gekrümmten Stoßlinien können dann durch Gerade 
ersetzt werden und die Zustände in den durch die 
Stoßlinien abgegrenzten Bereichen als konstant an- 
genommen werden. Wiein Bild1 dargestellt ist, gehen 
von dem Gabelpunkt die drei Stoßfronten A, a, b und 
die Unstetigkeitsfläche ce ab. Die Zustände in den vier 
hierdurch abgegrenzten Bereichen werden mit den In- 
dizes 1...4 gekennzeichnet, während die Indizes h, a, 
b sich auf den Hauptstoß und die beiden Teilstöße 
beziehen. : 

Ein einfacher Verdichtungsstoß ist vollständig 
bestimmt, wenn die Machsche Zahl der Anströmung M 
und der Stoßwinkel y oder zwei andere passend ge- 
wählte Parameter vorgegeben sind. HS feldes in der Umenvung dus Vers ER 
ist die Machsche Zahl der Anströmung M, für den eines gegabelten Verdichtungssto £ 
Hauptstoß und den ersten Teilstoß gleich. Ferner sind SR re ee re 
die Zustände nach dem ersten Teilstoß und damit M, und 4 ab. 
durch diesen bereits vollständig bestimmt, so daß IB 
beim gegabelten Verdichtungsstoß nur noch über vier Größen verfügt werden kann. An der 
Unstetigkeitsfläche C sind aber zwei Bedingungen zu erfüllen, nämlich Gleichheit des Druckes 
und der Strömungsrichtung, während Geschwindigkeitsbetrag, Dichte, Temperatur und Entropie 
unstetig sind. Insgesamt bleiben also zwei freie Parameter übrig, die den Gabelstoß bestim- 
men. Wenn einer der drei Stöße eindeutig vorgegeben ist, so ist damit der Verzweigungsstoß 
eindeutig bestimmt oder läßt wenigstens nur eine endliche Anzahl von Lösungen zu. 


a 


M2 ,P2 b 


2. Allgemeines Lösungsverfahren. 
Die Bedingungen des einfachen Verdichtungsstoßes werden bekanntlich dadurch gewon- 
nen, daß man auf die Stoßfront den-Energiesatz, Impulssatz und die Kontinuitätsgleichung an- 
wendet. Für unsere Zwecke erweist sich die von F. Schubert [4] eingeführte Parameter- 


Bild1. Schematische Darstellung des Strömungs- P 


TITTEN BER 


darstellung der Verdichtungstheorie als besonders vorteilhaft. Nach Schubert ist das Ver- 
hältnis des Druckes p’ nach dem Stoß zu dem Druck p vor dem Stoß gegeben durch: 

RL ER telze Ie  .. (1), 
wobei wir für die Klammer zur Abkürzung x gesetzt haben. Wenn x gegen 0 strebt, geht der. 
Stoß in eine Machsche Welle über. Auch das Verhältnis der Temperaturen, Dichten und die 
Entropie sind nur Funktionen von x. Dagegen ist der Ablenkungswinkel ö eine Funktion von 
M und «: 
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—— 


| s=arcctg | (SE me) Vers Ma)... 0: ex (2). 


M—(1+2) 


Aus der Theorie des einfachen Verdichtungsstoßes brauchen wir noch eine dritte Beziehung für 
die Machsche Zahl M’ nach dem Stoß, die hier abgeleitet werden soll, weil sie sich in der 
Arbeit von Schubert [4] nicht findet. Der Energiesatz lautet nämlich, wenn die Adiabate 
zugrunde gelegt wird: 
dp w: ee RE I ae 
3 a — — mo — t 
J3 Ba a An ri Ze un 


. Be w 
oder wenn wir die Schallgeschwindigkeit a? = x P_ einführen und ge M setzen: 


2 
EL ed 
= | a? T 
M? + a 
„—1 


Das Temperaturverhältnis 7/7’ ist aber bei Schubert angegeben: 


ET VE Hi! 
ee: 5: 


Damit haben wir alle Beziehungen aufgestellt, die wir aus der Theorie des einfachen Ver- 
dichtungsstoßes brauchen. Beim Gabelstoß treten als neue Bedingungen die Gleichheit des 
Druckes und der Strömungsrichtung an der Wirbelfläche c hinzu. 


PB = Pi 
oder 
Pa Pa 3 P4 (5) 
ee 
oder mit der Abkürzung A = log p’/p: 
= Act AR ee 
ON e u + N, Ta Re a re ee An» 


Nach Gl. (1) und (2) ist ö6=öd(M,x) =öd(M, A). Wir zeichnen die Schar der Kurven 
M = const in einem ö, A Koordinatensystem. Es ergeben sich hierbei herzförmige Kurven 
(Bild 2). Die praktische Durchführung dieser Auftragung wird dadurch wesentlich erleichtert, 
daß die Produktionswerte im wesentlichen bereits in der Arbeit von Schubert tabuliert 
sind (durch Verwendung von halblogarithmischem Papier kann man sich das Logarithmieren 
ersparen). Ferner benötigt man noch die Funktion: 
M?+ a a 
“zei Does a 
Re A 
u | 
die ebenfalls der Arbeit von Schubert entnommen werden kann. Damit ist man imstande, 
auf graphischem Wege Lösungen zu gewinnen. Wenn dieMachsche Zahl der Anströmung M, 
gegeben ist und der erste Teilstoß dadurch bestimmt ist, daß etwa M, festgelegt wird, findet 
man aus (8) A,. Man trägt ‘die Schar der herzförmigen Kurven zweimal auf, das eine Mal 
auf durchsichtiges Papier. In der einen Kurvenschar sucht man die Kurve M ı = const und be- 
stimmt den Punkt A, (Bild 2). In diesem Punkt trägt man die Kurve M, = const ab, indem 
man etwa das durchsichtige Kurvenblatt so verschiebt, daß der Nullpunkt mit diesem Punkt 
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_ zur Deckung kommt. Man erkennt aus Bild 2, daß der Schnittpunkt der Kurven M, = const 


und M, = const tatsächlich die gesuchte Gabellösung darstellt, denn dieser Schnittpunkt erfüllt 


‚die Bedingungen (6) und (7). 

Auch Eggink und Weise haben ein en 
graphisches Verfahren angewandt, das aber we- 
niger einfach als das hier entwickelte ist.” 

Es wird sich die Frage erheben, ob es außer 

dem trivialen Schnittpunkt A,, 6, stets einen 
zweiten oder gar mehrere Schnittpunkte der’'herz- 
förmigen Kurven gibt. Man wird finden, daß 
es je nach der Lage von A,, 6. keinen, einen oder 
zwei Schnittpunkte gibt. Die Bereichsgrenzen 
sind eng verknüpft mit den Grenzlagen des Ga- 
belstoßes, die wir im nächsten Abschnitt behan- 
deln wollen. Zuvor sei noch auf zwei Arten der 
Darstellung der Ergebnisse hingewiesen, die 
von Weise und Eggink vorgeschlagen 
worden sind. Das Ablenkungsdiagramm nach 
Weise, das vom Verfasser in verschiedenen 
Punkten berichtigt wurde, gibt eine sehr über- = 20° 10 0 10° 20°  d 
sichtliche Darstellung für die Ablenkungswinkel Bild. 2. Graphisches Lösungsverfahren zur Berechnung 
(Bild 3). Aufgetragen wird die Ablenkung du Tun ueniiönn, Aura It de Loknalhus der 
des ersten Teilstoßes über dem Ablenkungswin- beispiel: M, = 2,0, M, = 1,5. Die Gabellösung ist durch 
kel Ör, des Hauptstoßes. -Auf der 45°-Linie ist den Schnittpunkt der beiden herzförmigen Kurven 
6a = Ö1, d.h. 6,= 0. Diese Linie kann also einem Bern 
zweiten Teilstoß mit verschwindender Stoß- 
stärke (Machsche Welle) oder einem senkrechten zweiten Teilstoß zugeordnet sein, während 
die Waagerechte ö, = 0 einem ersten Teilstoß mit verschwindender Stoßstärke zuzuordnen ist 
und die Senkrechte ö, = 0 einen senkrechten Hauptstoß bedeutet. Links oberhalb der 45°- 
Linie liegt das Gebiet der gegendrehenden Gabeln, d.h. die Ablenkungsrichtungen der beiden 
Teilstöße sind hier entgegengesetzt gerichtet, während sie rechts unterhalb dieser Linie gleich- 
gerichtet sind. Im Ablenkungsdiagramm sind eine Reihe Kurven M, = const bzw. M,* = const 
eingetragen, wobei M,* mit der kritischen Schallgeschwindigkeit gebildet wird.!) Als typisches 
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ds 


i i i = ögli £ ie Ablenkungswinkel ö7, des Haupt- 

Bild 3. Übersicht aller bei Luft (% = 1,405) möglichen Gabeln, dargestellt durch die ‚bl ‚des Hz 

stoßes-und da des ersten Teilstoßes. Die Parameter der Kurven bedeuten die auf die kritische Schallgeschwindigkeit be- 

zogene Machsche Zahl der Anströmung M *, Im schraffierten Bereich existieren keine Gabellösungen (berichtigtes und 
ergänztes Diagramm nach A. Weise). 


‘spiel betrachten wir die Kurve M,*= 1,5. Sie beginnt an der ‚‚vorderen Grenzlage“ mit 
nn an 6, = — 11°, wobei der ‚erste Teilstoß in eine Machsche Welle übergeht. Die Kurve 
durchläuft dann den Bereich der ‚„Rückwärtsgabeln‘, bei denen der Hauptstoß stromaufwärts 
gerichtet ist, bis schließlich längs der Senkrechten DF der Hauptstoß senkrecht wird. Bei 
wachsendem ö, wächst auch ö, weiterhin an, bis schließlich längs der 45°-Linie DG der zweite 


(« +1) M3 
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2) M°:* = 


Zi 3 et a 3, wieder ab. Längs BG 
e- überschritten und ‚schließlich an der 45°-Li 


wird später noch ek sein. be: 
Von Eggink wurde die Auftragung im Bı usemanns 
 vorzugt (Bild 4). Aus dem Stoßpolarendiagramm kann die 
Hmacl nr und Richtung sowie auch der Stoßwinkel y ne: werden 

 auf.der Stoßpolare M,* = legt 
fest. Wir drehen nun M,* en 

Winkel ö, und zeichnen die Stoßpolare” M, 

Der. Zustand M,* nach dem zweiten zen m 

dieser Stoßpolaren liegen, ist aber nach 
Betrachtungen nicht mehr frei wählbar. I 
daher in das Stoßpolarendiagramm noch F 
zeichnet, die M,* = const. zugeordnet sind. 
punkt dieser M,*-Kurve mit der M,*-Kurve eı Mi. 208 
Den Punkt M,* findet man leicht dadurch, daß man PR Er 
BR ö, aus dem Diagramm abliest und den Winkel 
—=6,+ 6, vom Pol aus abträgt. Der Schnittpunkt B oe 
Bila-4 Die Aufregung im Busemann- Er Geraden mit der M,*- Stoßpolaren ergibt M,* (wo- 2 a 
schen Stoßpolarendiagramm (schematisch). bei zu beachten ist, daß nur der eine Schnittpunkt pRyERE 
kalische Realität besitzt). 2 >: 


3. Die ee des gegabelten Verdichtungsstoßes. . 
Das Ablenkungsdiagramm von Weise hat bereits erkennen lassen, daß der Hauptstoß 
bei konstanter Machscher Zahl der Anströmung in einem Bereich variieren kann, der dadurch 
begrenzt ist, daß auf der einen Seite der vordere, auf der anderen Seite der hintere Teilstoß in 
eine Machsche Welle übergeht. Bei abnehmender M achscher Zahl der Anströmung rücken 
beide Grenzlagen immer mehr zusammen, bis schließlich alle drei Stoßfronten nMachschen 
Wellen entarten. Als Existenzbereich der Lösungen kann man also im Stoßpolarendiagramm 


A 
ne 3 
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für den Hauptstoß ein krummliniges Dreieck erwarten, das begrenzt ist durch: i Rn 
Be -. 1. 2,=0 Vordere Grenzlage : 
Er, j 2. %,=(0 Hintere Grenzlage 
E- 3. M >» 
e* x ist dabei die in Gl. (1) gebrauchte Abkürzung. 
3 Er Dieses einfache Bild wird allerdings noch einige Abänderungen erfahren müssen. 
Be. .Zur Berechnung der Grenzlagen gehen wir von Gl. (5) aus, die wir unter Berücksichtigung 
% von a) schreiben: 
E.. | ee )- (Er tn: )- 
E Läßt man hierin x, bzw. x, gegen Null gehen, so erhält man nur das triviale Ergebnis 2, = x, 
Ber, bzw. &,= %,. Um weitergehende Aussagen zu gewinnen, muß man einen Grenzübergang vor-_ 
2 £ nehmen: : 
ni n—% 2% 
( 1 »_ 14 en > (9a) 
1 bzw. S 
d(, —%, 2x 
3 en —=1 + 771 EN a (ab), 


Für einen Sonderfall, nämlich den senkrechten Hauptstoß, läßt sich die Lösung für die 
vordere Grenzlage leicht angeben. Wegen cos y„ = cos y„—= 0 folgt nämlich in diesem Fall: 


dx d 
Fra (M? sin? Yyn—1) = 2 Mtsiny,cosy—=0 


dx, d 
— = —— 2 csın? Ber 
dx, EL LE 
Aus (3) und (4) ergibt sich aber durch Differenzieren nach x, an der Stelle x, — 0): 
dx, =0) 1 = A| rk Vo Ar (10). 


ED 


er ee orrel ist in etwas ee Form bereits von A. We ise[5] angegeben neh, 


& stoß M,=1,484. - 
- Für den allgemeinen Fall der Grenzlagen sind die Gl. (98) und (9b) nicht ausreichend, 
Wir müssen auch Gl. (7) heranziehen, die wir in der Form schreiben: 


£ E = | n x ör ( (M;» 2) Ön (M;, &,) sE ö,(M;, %;) Br .. # . (12). 


Um zu den Gleichungen der beiden Grenzlagen zu gelangen, müssen wir auch hack nach 
0. % bzw..x, differenzieren, wobei im einen Fall wiederum «, = 0, x, = x, und im anderen Fall 
0, % = 2, zu setzen ist. Man erhält damit folgende beiden Gleichungen für die beiden 
..... Grenzlagen: | 


ee; | dd, da, (43% | (22) da, (2% ) am: 

Br i da, de, a oo,  \0% ee dx, oM?) dx, 

4 ya : a ER (H=) . 

= x ; - R P : (M,=Mı) - 

z ee ee ee 

m da, da, 2 (2, ) "da, en oM? 1 3 I 
(2, =%p) 6:7: 
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 “ Nach Ausführung der Differentiationen und Einsetzen des ee erhält man für die 
vordere er die Gleichung: 


2 2. ym—I if 2 ) | 
Plt+ 47° je ee IB = 


und für die hintere Grenzlage: 
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Aus den Gleichungen sa und (14) kann zu jeder Mach ones Zahl M, der Anströmung 
ein bestimmtes x, = M? sin 2y, — 1, also auch der Stoßwinkel y,„ berechnet “werden, bei dem 
gerade einer der beiden Teilstöße in eine-Machsche Welle übergeht. Für eine gewisse Mach- 
sche Zahl entarten alle drei Stoßfronten in Machsche Wellen. Um dieses M, zu gewinnen, 
müssen wir in (13) und (14) x,= 0 setzen. Es zeigt sich aber, daß man nur über einen 
weiteren Grenzübergang diese Machsche Zahl gewinnen kann. Schreiben wir nämlich (13) 
in der Form F (M,, x,) = 0 und entwickeln in eine Ta ylorsche Reihe an der Stelle x, = 0: 

2 
Seerde Bes 0 A SS, 
% % 


so sind die beiden ersten Glieder Null, und erst das Verschwinden’ der zweiten Ableitung liefert 
die gesuchte Bedingungsgleichung: 


RR F(M,=)=F(M,0) + 


‚M5—3%) HM? (12 44) +8=0 ........ (15), 
welche für x = 1,405 die beiden Lösungen: 
M, = 1,2456 


und M, = 2,5648 liefert. 


er Sie liefert für Luft (x = 1,405) als untere Grenze eines Gabelstoßes mit senkrechtem eg “= 


NR 


g RT 3 4 > 1 2 
Man erken die ı ne re u 


längs der Kurve ABDE der vordere, längs A@ der hintere Teilstoß : 
über. Bei y„ = 90° wird der Hauptstoß senkrecht, was der Linie DF 
„Oberhalb DF liegt das Gebiet der Rückwärtsgabeln, wo der Hauptst 


a 


Bild 5. Die Grenzlagen des Gabelstoßes für einen gegebenen Hauptstoß betrachtet. Im schraffierten Bereich gibt es zu einem { 
g>z2banen Hauptstoß zwei, im kreuzschraffierten Bereich drei Gabellösungen. Oberhalb der strichpunktierten Linie führt . a 
der Hauptstoß auf Unterschallgeschwindigkeit. s . . 


unterhalb DF das Gebiet der Vorwärtsgabeln, wo der Hauptstoß stromab weist. Unterhalb der . 
strichpunktiert eingezeichneten Linie BG führt der Hauptstoß auf Überschallgeschwindigkeit 
(M,> 1). Beim Punkt A schrumpft der Lösungsbereich zu einem Punkt zusammen. Alle drei he 
Stoßfronten entarten in Machsche Wellen. Dieser Punkt entspricht der unteren Gabel- i 
grenze, d.h. unterhalb einer Machschen Zahl der Anströmung M, = 1,2456 (für «= 1,405) 
ee -. sind Gabelstöße mit den Stoßbedingungen nicht vereinbar. Ein wesentlich anderes Bild bietet 
e: der zweite Lösungsbereich, der erst bei Machschen Zahlen M, > 2,5648 (für x = 1,405) ver- 
= wirklicht werden kann. Variiert man nämlich y, bei konstantem M, von der hinteren Grenzlage - 
ausgehend, so erreicht y, ein Maximum, ohne daß die vordere Grenzlage erreicht ist, der man 
sich erst bei abnehmendem y, nähert. Im schraffierten Gebiet gibt es daher zu einem gegebenen 
Hauptstoß zwei Gabelstöße, bei Berücksichtigung des ersten Lösungsbereiches sogar im kreuz- 


= schraffierten Gebiet drei Lösungen. 
3 Um den Vergleich mit dem Ablenkungsdiagramm (Bild 3) zu erleichtern, sind die charak- 

4 teristischen Punkte mit gleichlautenden Buchstaben bezeichnet. Der zweite Lösungsbereich ist 

y 

E 


Bild. 6. Auftragung der Grenzlagen im Stoßpolarendiagramm für den Hauptstoß (links) und den zweiten Teilstoß (rechts). 

Die Buchstaben entsprechen denen des Bildes 3 und 5. Faßt man das linke Diagramm als Stoßpolare des ersten Teilstoßes 

auf, so muß dieser stets auf Punkte rechts von AG führen, während rechts der Linie Uy zu einem gegebenen ersten Teil- 
stoß zwei Gabellösungen gehören. 
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hier auf den Winkelraum der gleichdrehenden Gabeln beschränkt. In Bild 6 sind die Lösungs- 
bereiche auf das Stoßpolarendiagramm für den Hauptstoß und den zweiten Teilstoß über- 
tragen. Die Grenzkurve der unendlich großen Machschen Zahl der Anströmung wird hier. 


durch die Kreisberandung gebildet. 


4. Vergleich mit Messungen. 


Gegabelte Verdichtungsstöße sind wohl zuerst bei Lavaldüsen beobachtet worden, bei denen 
der Gegendruck so groß ist, daß die Strömungsvorgänge nur durch Verdichtungsstöße realisiert 
werden können. Verdichtungsstöße, die von Überschallgeschwindigkeit auf Unterschallgeschwin- 
digkeit führen, lösen aber erfahrungsgemäß beim Auftreffen auf eine feste Wand die Strömung 
ab. Von der Ablösestelle an der Wand geht offenbar ein schräger Verdichtungsstoß ab, da nur 
dieser eine Ablenkung der Überschallströmung von der Wand weg verursachen kann. Wegen 
der Kleinheit dieser Ablenkung geht der Schrägstoß nahezu unter dem M ach schen Winkel ab. 
In seinem weiteren Verlauf kann der Schrägstoß zwei wesentlich verschiedene Formen an- 
nehmen: 

1. Der in Wandnähe beginnende Schrägstoß geht mit zunehmendem Abstand von der Wand 
in ein senkrechten Stoß über. Diese Stoßform verursacht nur geringe Ablenkung und 
zeigt eine starke Neigung zum turbulenten Wiederanlegen der Strömung. Eggink[6]hat 

.diesen Fall bei M, = 1,2 und 1,3 gefunden. 

2. Schlierenbeobachtungen von Weise[5] und Eggink[6] haben jedoch gezeigt, daß bei 
größeren Machschen Zahlen eine Verzweigung des Stoßes eintritt. Weise[5] hat ober- 
halb der theoretischen Gabelgrenze (M, = 1,2456) nie Verdichtungsstöße ohne Verzweigung 
gefunden. Den engen Zusammenhang mit der Ablösung der Grenzschicht hat er dadurch 
erwiesen, daß durch Grenzschichtabsaugung die Gabelbildung verhindert werden konnte. 

Eine neuere experimentelle Arbeit aus dem Institut von Prof. Ackeret[7] scheint im 
Widerspruch hierzu zu stehen. In dieser Arbeit wird die Vermutung ausgesprochen, daß bei ' 
laminarer Grenzschicht ein gegabelter Verdichtungsstoß auftritt und zwar bei Machschen 
Zahlen, die wenig über 1 liegen, Vielfachstöße, die bei Erhöhung der Geschwindigkeit an Zahl 
abnehmen, bis schließlich ein einfacher ‚‚A-Stoß‘ erreicht wird. Dagegen ergeben turbulente 
Grenzschichten nach Ackeret nur stetig gekrümmte Verdichtungsstöße ohne Verzweigung. 
Dabei haben aber Ackeret und seine Mitarbeiter in einem Bereich der Machschen Zahlen 


. experimentiert, in dem Gabelstöße mit den Stoßgesetzen nicht vereinbar sind. Obwohl die bei 


turbulenter Grenzschicht beobachteten Doppelstöße in der Züricher Arbeit mit Kondensations- 
erscheinungen in Zusammenhang gebracht werden, erscheint eine solche Deutung für die beob- 
achteten Stoßvorgänge wenig wahrscheinlich, weil der Züricher Kanal mit einem Gebläse be- 
trieben wird, die Luft also trocken ist.)) Die Druckverteilungsmessungen lassen auch die 
Deutung zu, daß der auf die Wand zulaufenden Stoßfront ein keilförmiger Verdichtungsfächer 
vorausgeht, der auf den Schlierenbildern einen Gabelstoß vortäuscht. 

A.FageundR.F, Sargent[8] untersuchen in einer experimentellen Arbeit den Ein- 
fluß einer turbulenten Grenzschicht an einer geraden glatten Wand auf senkrecht und schräg 
auffallende Verdichtungsstöße. Der Überschallkanal ist dabei so ausgebildet, daß nur die Gegen- 
wand nach Art einer Lavaldüse gekrümmt ist. Die annähernd senkrecht auf die Wand auf- 
treffende Stoßwelle zeigt zwei Erscheinungsformen. Während die Stoßwelle bei den Mach- 
schen Zahlen M, = 1,08 und 1,206 ohne Verzweigung in die Grenzschicht hineinläuft, beob- 
achtet man bei M, = 1,310 am Rande der Grenzschicht eine schwache Andeutung einer Gabel- 
bildung. Mit wachsender Machscher Zahl der Anströmung rückt der Gabelpunkt zunächst 
schnell, dann immer langsamer in das Innere des Strömungsfeldes und erreicht bei M, = 1,443 
eine Entfernung von 12 mm von der Wand. Gleichzeitig wächst dabei das Druckverhältnis 
?,/p, an, und P,/P, = 1,8 scheint die Rolle eines kritischen Druckverhältnisses zu spielen. Diese 
Beobachtungen bestätigen also die Auffassung, daß Gabelstöße erst von einer gewissen Mach- 
schen Zahl der Anströmung an existieren können. 

Durch Einfügen eines Keiles an der Gegenwand wird erreicht, daß von dort ein schräger 
Verdichtungsstoß ausgeht und an der geraden ‚‚Meßwand‘“ reflektiert wird. Sobald jedoch der 
Einfallswinkel einen kritischen Wert überschreitet, löst sich die Stelle, wo die Reflexion statt- 
findet, von der Wand ab und rückt in das Innere des Strömungsfeldes. Dabei bildet sich eine 
Art „umgekehrter Gabelstoß‘“, wobei die beiden Teilstöße durch den einfallenden und reflek- 
tierten Stoß gebildet werden, während eine dritte Stoßwelle zur Wand läuft. Eine solche Stoß- 
form ist erst bei genügend großer Machscher Zahl der Anströmung möglich. Durch Änderung 
des Einfallswinkels kann die Stoßstärke geändert werden. Während bei p,/P, = 1,76 und 


»)L.Prandtl hat die Vermutung ausgesprochen, daß es sich bei den Vielfachstößen vielleicht um eine gekoppelte 
Sohwingung des Systems Druckfeld — Grenzschicht handelt. 


" glückte Versuch einer Theorie des Gabelstoßes unternommen. Da hierbei im Üh 
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[1] iR Weise: Theorie des gegabelten Verdichtungsstoßes. Jb. deutsche Luftfahrtforsch. 1943 IT AO 


‘ Ringegangen am 23. 7. 1947. 
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M, = 1,439 die oben beschriebene Stoßlorm beobachtet worden ist, gal 
und M, — 1,439 die zur Wand laufende dritte Stoßwelle nochmals nach 4 
die Wand auftreffenden Stoßwellen, so daß ein doppelter Gabelstoß entsteht 
gungspunkten. Übrigens wird von Fage und Sargent auch der allerdi 


keine untere Grenze der Machschen Zahl für die Existenz von Gabelstößen in E 
tritt, muß das Fehlen von verzweigten Stößen bei kleinen M a chschen Zahlen lure x 
schichteinflüsse plausibel gemacht werden. . ME e ee a; 
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Ein Beitrag zur Theorie der Wellen in freien Gasstrahlen. 
Von Fritz Kai in Wittenberg. 


Es wird nachgewiesen, daß im Innern ebener Gasstrahlen schiefe Verdichtungsstöße auftreten können, 
die die Symmetrie des Strömungsvorgangs stören. Daraus wird gefolgert, daß die bisherige Auffassung, 
die eine Ausbildung regelmäßiger Wellen in solchen Gasstrahlen annimmt, nicht allgemeingültig sein kann. 

It is demonstrated that inclined shocks of compression may appear in the interior of plane gas jets, and 
disturb the symmetry of the current. Hence it is concluded that the opinion hitherto existing cannot prove 
right in general, namely that regular waves are arising in such gas jets. : 


On demontre qu’&a l’interieur d’un jet gazeux peuvent prendre naissance des chocs de compression inclin&s 
qui detruisent la symmetrie du phenomene d’ecoulement. On en deduit que l’opinion admise jusqu’ä present, 
qui admet l’existance d’ondes regulieres dans les jeus gazeux de ce genre, ne saurait &tre applicable d’une 
fagon generale. 

lokassIBaeTca, YTO BHYTPH INIOCKHX TASOBEHX CTPyÜ MOoryT HaÖMONATECHA HAKNIOHHEIE 
CIKATHA YAaPHOTO XapaKkTepa, KOTOPHe HAPYIHAIMT CHMMETPHIO ABJIeHHAd TeyeHun. OTcioya 
BEIBOAHUTCH, YTO TOCIHONCTBOBABIIeEE MO CHX TOP MHEHHe O HAJIHYHH PABHOMEePHEX BOAH B 
HONOÖHEIX TA30BEIX CTPYAX He MO3KeT ÖBITB CHPABE/KIHBEIM BO BCEeX CIIyyasıx. 


Seit langem sind die Wellen experimentell!) bekannt, die in dem Strahl auftreten, der sich 
beim Ausströmen eines Gases mit Überschallgeschwindigkeit aus einer Mündung bildet. Mit den 
Methoden der Gasdynamik ist die Strömung eines solchen Strahls für einen ebenen Strömungs- 
vorgang auch theoretisch untersucht worden, worüber L.Prandt] zusammenfassend in seiner 
„Strömungslehre‘“) berichtet. Dort wird die Auffassung vertreten, daß in einem Gasstrahl, der 
als ebene adiabatische Potentialströmung angesehen wird, sich regelmäßige Wellen ausbilden. 
Der Begriff ‚regelmäßige Welle“ wird nicht näher definiert; doch ist er wohl in der Weise zu 
verstehen, daß die Strahlglocken, wie ich die entstehende Strömungsfigur (vgl. Bild 2) bezeichnen 
möchte, geometrisch und energetisch völlig symmetrisch sind. Begründet wird diese Auffassung 
damit, daß am Ende der ersten Strahlglocke derselbe geometrische und energetische Zustand -. 
herrscht wie anfangs in der Mündung (vgl. Fig. 233 l.c., wo der Druck an beiden Stellen mit 
p, angegeben ist) und dadurch die Voraussetzung für eine fortdauernde immer gleiche Wieder- 
holung des Strömungsvorganges geschaffen wird. . 

Bei gasdynamischen Arbeiten, die am Lehrstuhl für technische Mechanik der Technischen 
Hochschule Dresden unter Leitung von W. Tollmien durchgeführt wurden, ergaben sich 
Zweifel an der allgemeinen Gültigkeit dieser Auffassung. Jedoch konnte damals aus zeitbedingten 
Gründen dieses Problem nicht eingehend untersucht werden. Eine überschlagsmäßige Untersuchung 
schien allerdings das von Prandt langegebene Ergebnis zu bestätigen ; doch konnte eine bündige 
Entscheidung nicht getroffen werden. 

Mit den in Dresden entwickelten methodischen Hilfsmitteln habe ich nun das Problem der 
Wellenbildung von freien Gasstrahlen erneut angegriffen und bin zu einem teilweise von der 
obigen Auffassung abweichenden Ergebnis gelangt. 


fi 
N. 2 L 


‘) E. Mach und P. Salcher, Sitzungsber. d. Wiener Akad., math.-naturw. Klasse 98 IIa (1889), S. 1303. 
’) L. Prandtl, Führer durch die Strömungslehre, 3. Aufl. Braunschweig 1942. S. 255-257. 


Bei der Konstruktion des Strömungsfeldes ergab sich nämlich im Innern der Strahlglocke 
. ein Zusammenlaufen gleichartiger Ma chscher Wellen, die dabei augenscheinlich eine Enveloppe 
bildeten. Eine solche Enveloppe stellt aber nach Untersuchungen von W. Tollmie n®) eine 
 Grenzlinie der adiabatischen Potentialströmung dar, über. die hinaus eine Fortsetzung der an-- 
kommenden Strömung nicht möglich ist. Vielmehr bricht .die Potentialströmung bereits strom- 
aufwärts zusammen, und an der Stelle entsteht in der Strömung ein mit verschwindender Stärke 
er schiefer Verdichtungsstoß, durch den offenbar die bis dahin bestehende Symmetrie 
‚gestört wird. 


Zur Konstruktion derartiger Überschallströmungen, die schiefe Verdichtungsstöße nieht- 
konstanter Anströmung enthalten, ist ein neues auf differentieller Grundlage beruhendes Ver- 
fahren entwickelt worden. Zu dem Ende werden der Ablenkungswinkel ö der Stromlinie sowie 


der Geschwindigkeitsbetrag w' hinter dem Stoß als Funktionen des Stoßfrontwinkels y und des 


Geschwindigkeitsbetrages w vor dem Stoß aufgefaßt 4) ä 


ER 2 | N a a er (1). 
Bildet man die totalen Differentiale von ö und w 
06 06 ow dw’ ; 
dö zz — —— je x sen — 
FB GR 5, dw; dw 37 dy-+ a EEE a (2) 
und eliminiert 'dy aus diesen Gleichungen, so erhält man, wenn man noch durch . 
86 . 10w 

A 0y H B= up ee et Snake ler tie Alter MiRh esse (3) 

und 
CR) 1 ow' = 
( FR a en a ee (4) 


die Hilfisgrößen A, B, C, D einführt, die im folgenden als ‚‚Stoßfrontgleichung“ bezeichnete 
Beziehung 

vn A dw‘ (A D 

- Ba z +oldu=0 BT RE RE ER (5), 


welche die im allgemeinen unbekannten differentiellen Änderungen des Ablenkungswinkels ö 
-und des Geschwindigkeitsbetrages w’ hinter dem Stoß mit der bekannten Änderung des Geschwin- 
digkeitsbetrages w vor dem Stoß längs der Stoßiront miteinander verknüpft. Aus Gl. (2) folgt 
für die zugeordnete differentielle Anderung des Stoßfrontwinkels y 
dö CC | 
re ep a DE Er N EEE Da Ba Pa (6). 
Die an sich hinter einer Stoßfront auftretende Wirbelbildung kann vernachlässigt werden, solange 
sich der Stoßfrontwinkel y nur wenig längs der Stoßfront ändert, d.h. also solange ihre Krümmung 
day 
do 
ten Art fast immer der Fall. Daher kann auch die Strömung hinter dem Verdichtungsstoß wieder 
als Potentialströmung angesehen werden, allerdings mit einer anderen Bernoullischen Kon- 
stanten als vor dem Stoß, was beiden Druckzahlen durch Berücksichtigung des Drosselfaktors K°) 
zum Ausdruck kommt. 


klein bleibt. (da Bogenelement der Stoßfront.) Das ist bei Strömungen der hier behandel- 


Um die Konstruktion des Strömungsfeldes hinter dem Verdichtungsstoß weiterführen zu 


können, muß die Stoßfrontgleichung (5) mit den vom Prandtl-Busemann -Verfahren 
her gut bekannten ‚‚charakteristischen Hodographen“ 


aan =0 RE ee (7) 


kombiniert werden. (9 Stromlinienwinkel, aM achscher Winkel) Ersetzt man in erster Nähe- 


rung, wie es bei den gasdynamischen Näherungsmethoden üblich ist, die Differentiale durch _ 


Differenzen, wobei man bei der Bildung der Differenzen längs der Stoßfront bzw. einer Mach- 
schen Welle von einem Punkt P, mit den bekannten Strömungsgrößen y. wg U, bzw. 
Il (wyr, dır) zu dem Punkt P mit dei unbekannten Strömungsgrößen y, w', 0° (vgl. die schemati- 


») W. Tollmien, Grenzlinien adiabatischer Potentialströmungen. Z. angew. Math. Mech. Bd. 21 (1941), S. 140 
bis 152. * € 

4) Die Anregung zu dieser Entwicklung verdanke ich Herrn M. Schäfer,; dem an dieser Stelle herzlich ge» 
dankt sei. ! 

s) Vgl. hierzu etwa F.Schubert, Zur Theorie des stationären Verdichtungsstoßes. Z. angew. Math. Mech. Bd. 23 
(1943), S. 129—138. 
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= er 


sierte Darstellung in Bild 1) fortschreitet, so Be sich die unbekannten Be 23 = 
aus den Gl. (5) und (7) bestimmen. Man erhält bei einer linksläufigen Stoßfront (y > 0) air Unis 5: 
läufigen Machschen Wellen mit der Neigung 9 +. für den Geschwindigkeitsbetrag w’ : ; 


+9 dt Ge etgon + IrE, a 5 
ee Ver 


4, „ Agaı 
Bow, wır 
und für den Stromlinienwinkel ee 


= dr — ctgarı+ = 


st Stoßfrent 
/MW linksi-Machsche Welle 


Bein Dero) 2 -40R, 


Bild 1. Die in den ee ee: vorkommenden ungestrichenen 

Strömungsgrößen sind als bekannt anzusehen, da sie ja mittels der be- 

kannten Näherungsverfahren zur Konstruktion adiabatischer ebener Potentialströmungen gewon- 

nen werden können. Der Stoßfrontwinkely im Punkt ?P wird im Rahmen der Näherung be- 
stimmt aus 


4 


6—ö6, , (0 
y=Yyo RER ra, ee a a er a A (10). 


Die Güte der so erhaltenen Näherungslösung kann in jedem einzelnen Falle mit der aus der Ver- 
dichtungsstoßtheorie bekannten exakt geltenden Gleichung?) 


x—1 1 
tg (y— ()— Ri rt) en (11) 


geprüft werden. (x= c,/c, Verhältnis der spez. Wärmen.) Dabei hat die praktische Erfahrung 
immer eine recht befriedigende Übereinstimmung ergeben. Auf diese Weise läßt sich die Stoßfront 
ebenso bequem punktweise aufbauen wie das übrige Strömungsfeld. 

Die Voraussetzung für die praktische Anwendung des beschriebenen Verfahrens bildet 
natürlich die Kenntnis der vier Hilfsgrößen A, B, ©, D. Diese lassen sich jedoch durch einfache 
Differentiationen aus den Hauptgleichungen der Verdichtungsstoßtheorie 


“+1 I | 
eh '&y | Marsny —T sin?y—1 x d) 
(Ma Machsche Zahl der Anströmung unmittelbar vor dem Stoß) 
und 
a DR 
er ee (ID) 


gewinnen. Nach Ausführung der Differentiationen [vgl. die Definitionsgleichungen (3) und (4)] 
ergibt sich nach einiger Rechnung: 
(#*+ 1) Ma* sinty— (Ma?sin?y — 1)? — 2 ı . Ma? (Ma?sin®y +1) 


Az — ———— 
a eg ae 
(Ma?sin®y — 1) (x. Ma?sin?2y +1) — = Matt/ a Ma?sin?y 


B=tgy—tgy—d)(l—A). Scare, See UT er (13), 
1 3 
aa (x — 1) Ma2 costy 87 
Om ———— Ta I ee a re es 
(Ma?sin®y —1) (x.Ma?sin?y +22 Ma ’ “EL Mas sin? 
1 
DSmt 90. Sa Re Sala) 


*) Vgl. hierzu und zu den folgenden Gl. (I) u. (II) aus der Verdichtungsstoßtheorie F. Schubertl.e. 


ces) ut, 


beträgt 


PR a2 


strahlen 


Rn > VEREIN a er 
r it EHOUNIEN rn 1 1en Gas 
” 2 v * Pr gr = er - - v2 


Nach diesen Ausführungen über die benutzten methodischen Hilfsmittel sollen nun die 


an den Beispielen gewonnenen Ergebnisse dargelegt werden. 


£ Ebener Gasstrahl Ma-150 Pulbr2 = 
Dh ‚Srömingfd 


= 


Sz2 
II Sa 
Bild 2. r 


Bei den beiden konstruierten Strömungsbeispielen (s. Bild2 und 4) handelt es sich um den 
'ebenen Strahl eines Gases (x = 1,405, etwa trockene Luft), der beim Ausfluß aus einer länglich 
rechteckigen Mündung gegen Unterdruck entsteht. Die Geschwindigkeit der Gasströmung im 


Mündungsquerschnitt entspricht in beiden Fällen der Machschen Zahl Ma= 1,50. Das Ver- 


hältnis vom Mündungsdruck ?y zum konstanten Außendruck p, der ruhenden Umgebung 


er PulPpa —=2,0 bzw. Pulpa —5,0. 


Die Bilder 2 und 4 zeigen die Strömungsfelder der Strahlen mit den eingezeichneten Mach schen 


Wellen. Die im Innern der Strahlglocken entstehenden Verdichtungsstöße sind durch stärkere 
Strichart hervorgehoben. In der gleichen Strichstärke sind auch die Freistrahlgrenzen gezeichnet. 
Im Bild 4 sind außerdem zur. Veranschaulichung des Strömungsvorgangs noch einige Stromlinien 
eingetragen, zwischen denen jeweils der gleiche Durchfluß besteht. 

Beim ersten Beispiel (Bild 2) sind zur Klärung des Strömungsvorgangs zwei vollständige 


Strahlglocken konstruiert worden, an deren Ende jedesmal ein mit verschwindender Stärke 


beginnender schiefer Verdichtungsstoß einsetzt. Durch diese Verdichtungsstöße wird die bis 
dahin vorhandene Symmetrie der Strömung gestört, was am besten ein Vergleich der geometrischen 
und energetischen Strömungsgrößen in der Mündung (Index M), am Ende der ersten (Index #1) 
und am Ende der zweiten Strahlglocke (Index E 2) zeigt. Die Breite des Strahls in diesen engsten 


- Querschnitten bezogen auf die Mündungsbreite 2b beträgt: 


bu/b=1,0; byb—=1,18; by,/b = 1,36. 


Die Öffnung des Strahls bestimmt durch den Abgangswinkel © der Freistrahlgrenze verläuft 
wie folgt: i 


Oel Heli, 9, = 12°32), 
und der mittels des Kesseldruckes p, dimensionslos gemachte Druck zeigt folgenden Gang: 


PulPo = 0,2722; PEıl/Po —= 0,2711; PralPo — 0,2696. 


Druck-und beschwindigkeitsverteilung langs Strahlachse 
Ma= 750; D,/D4°2 b 


Bild 3. 
6* 


Ü ‚Die Ausdrücke sehen zwar in der allgemeinen Form recht kompliziert aus, lassen sich aber für vor- 
Kr Et Werte von Ma und y in jedem einzelnen Falle ohne großen Aufwand numerisch an- 
ER 32) Te Fr ; | | 


5 


ig a  ı 
4 Kai, Bin Beitrag 
Besonders augenfällig kommt die g 


indi 


Beim zweiten Beispiel (Bild 4. Es ist nur die obere Hälfte des Strömungsfeldes gezeichnet, j 
die untere ist symmetrisch zur Achse zu denken) genügte bereits die Konstruktion einer ng u 
glocke zum Nachweis der Unsymmetrie der Strömung. Hier laufen nämlich die ersten an 

Freistrahlgrenze reflektierten M a c h schen Wellen bereits in der Mitte der Strahlglocke zusammen. 
Der hier sich ausbildende Verdichtungsstoß wird seinerseits an der Mittelachse, reflektiert und 
durchläuft, wobei er sich zu beträchtlicher Stärke aufschaukelt, die ganze Breite der Strahlglocke. 
Dadurch ist im hinteren Teil der Glocke die Strömung völlig verändert. Dies wird besonders 
deutlich an Hand von Bild 5, das wieder Druck- und Geschwindigkeitsverteilung längs der Mittel- 


Druck-und Geschwindigkeitsverterlung längs Strahlachse 
" Ma« 150; Py/pa"5 
100 


i 
| 


achse des Strahls darstellt. Man erkennt dort den infolge des Verdichtungsstoßes bei x/b = 23,3 
auftretenden Druck- und Geschwindigkeitssprung und sieht ferner, daß von dieser Stelle ab die 
Strömung von einem anderen Druck- und Geschwindigkeitsverteilungsgesetz beherrscht wird. 
Ein Vergleich der geometrischen und energetischen Zustandsgrößen in der Mündung mit denen 
am Ende der Glocke liefert einen weiteren Beweis für die Umgestaltung des Strömungsvorgangs. 
Man hat nämlich: 
und 
On = 28° 10°, GO; = 20° 25. 

sowie j 

PulPo = 0,2722; PR, = 0,1455. 


Daraus erhellt, daß die zweite Strahlglocke wegen des für sie völlig veränderten Ausgangszustandes 
keinesfalls symmetrisch zu der ersten werden kann. Auch hier ist also die Ausbildung einer regel- 
mäßigen Welle nicht möglich. 

Mit diesen vorstehenden Ergebnissen ist der Nachweis erbracht, daß die eingangs angeführte 
Auffassung einer Ausbildung von regelmäßigen Wellen in ebenen Gasstrahlen zumindest keine 


esit zwei konkreten ee ee ein im ae 77 
) eh eher werden. Hierdurch weicht der geometrische und 
ustand ‚am Ende einer Strahlglocke teilweise erheblich von,dem Ausgangszzustand 
r Mündung ab, und damit entfällt: die Voraüsseizung für eine fortdauernde immer 


; Wiederholung des Strömungsvorgangs. 


Da bekanntlich im Verdichtungsstoß Sirömungsenergie verzehrt us stellen die vorliegen- E 
en Strömungsbeispiele gedämpfte Wellen dar, die — auch ohne Berücksichtigung irgendwelcher 


gseinflüsse — nach einer endlichen Zahl von Strahlglocken abklingen und zur Ruhe kommen. 
Ob und wi wie weit überhaupt eine regelmäßige ad infinitum fortdauernde Welle möglich ist, bleibt 
ne noch offene Frage. Diese glaube ich Re nach meinen bisherigen Be a ae in negativem 
Sinne beantworten zu können. 
i ange: I Juni * 


Be tkeihoodSchätzine und Recke 
en der St atistischen Abteilung des W. G. Kerckhof f-Instituts.) 
> - Von Maria Pia Geppert in Bad-Nauheim. 
Es wird gezeigt, daß beim Rückschluß von einem Stichprobenwert auf einen Parameter im Kollektiv 


weder Dichtemittel noch Erwartungs- noch Zentralwert der Mutungsverteilung mit der „besten“ Schätzung j 


nach dem -Maximum-likelihood- Prinzip zusammenfallen müssen. 


Itis shown that, in estimating a parameter in a population from a statistic in a ‚sample, neither the mode, 
nor the mean value, nor the median of the fiducial distribution must coincide with the „best“ estimation 
aecording to ihe principle of mazimum-likelihood. 


. On demontre que, en estimant la valeur d’un parametre dans un collectif par un Echantillon, ni la nie: 
ni la moyenne, ni la mediane de la distribution fiduciale doivent coincider avec c la ymeilleure« estimation 
selon le principe de mazximum-Llikelihood. 

lokassiBaetcH, YTO IPH oMeHke mapamMmeTpa KOMIeKTUBA IIYTeM BSATUA OTAENBHHX IPOÖ 
Hu CpeAHAA INIOTHOCTb HH O3KHNACMOR HIN IEHTPA.IBHOe SHAYEHHE BEPOATHOTO PacHperelteHun 
He AOIBEHEL o6A3aTeJIbBHO COBHANATB C ‚„‚„HauIyymeii‘‘oNeHkof Io MeTony ‚maximum likelihood‘, 


. Sei © ein Parameter einer Gesamtheit (z.B. Medianwert, Streuung oder Korrelations- 
koeffizient derselben usw. ), 7 der entsprechende kontinuierlich veränderliche Parameter in Stich- 
proben bestimmten Umfangs, der zur Schätzung des unbekannten Gesamtheitsparameters © 
dient. Die kontinuierliche Wahrscheinlichkeitsverteilung von Tin Stichproben gleichen Umfangs 
bei festem © sei gegeben und ihre Dichte laute 

/(T|0) =0 a BL ae TE 3 
= Das — mit gewissem Recht meist nach Bayes benannte — Fückschlufproblem der 
mathematischen Statistik (engl. ‚problem of estimation‘“) lautet: „Was läßt sich auf Grund 
_ eines in einer Stichprobe beobachteten T-Wertes über den unbekannten Wert des zugrunde 
liegenden © aussagen?“ Zu seiner Lösung sind bekanntlich im Laufe des letzten halben Jahr- 
hunderts — zum Teil aus wesentlich älteren Wurzeln — vier Methoden entstanden!), die zwar 
zum Teil eng ineinander verflochten sind, aber doch auf verschiedene Fragen Antwort zu erteilen 

suchen (vgl. [2a]) und hier unter Verzicht auf nähere Begründung kurz genannt seien. 
I. Nach dem Prinzip der Maximum-likelihood wird als ‚beste‘ Schätzung von © derjenige 

— von T abhängige — Wert 


FEB a ET ey 
Yehzler. für den die Wahrscheinlichkeitsdichte f(7|©) maximal ist, für FAR al 
of(T|®) 
FT) Or nei) 
e=LT) 
f(T|O 
mn NEN WEN RE (4,) 
0-7) 
oder allgemeiner mit einem bestimmten k>1 
”f(710) ee aper..n.n, 0.0) 
20 
e=L(M) 
#*f(T]0) 
ae < 0 ee een Be arte te Ve I (4;) 
8=L(T) 


gilt. 
1) vgl. die Zusammenstellung in [2a] und die dort zitierte Literatur. 


om zwei ERREN Be ee Be 
ge mit I Aussage SER 


von © aus nach unten bzw. oben abgegrenzte äufallsbepeiche 2” F- gerade og den beo} bach! 


- Mutungswahrscheinlichkeil iv (a T) d© mit der Dichte En SE 


$- EN Ihre genaue eos ist in [2b], S. 10, auseinandergesetzt. Den deckt sich a 


er Yin) = 9<029- 


zu irren, bezogen auf alle möglichen T- und O-Werte, 2E Bas Diese Mut u 
demnach die extremsten @-Werte, deren der Gesamtwahrscheinlichkeit (1-2) ents 


T-Wert enthalten. Ir Be 


$ 


"III. Nach R. A. Fisher [1a, b] läßt sich Drama dlas kontinuierliche en 


w(lT) = frel® a. ee ITS ee % 
nu 


pi 


die ’Angabe dieser Rückschlußwahrscheinlichkeit inhaltlich genau mit der Angabe der Mutungs- K 
grenzen Gl. (6) als Funktion der veränderlich gedachten Zufallsziffer und ist letzterer — allerdings ra 
auf Kosten der Anschaulichkeit — nur formal überlegen. 3 r er 


IV. Verändert man die Fragestellung von Grund aus, indem man sich die Gesamtheit F 
mit dem unbekannten Parameter © zufallsmäßig einer Obergesamtheit solcher Gesamtheiten 
entnommen denkt, innerhalb welcher @ nach dem a-priorischen Wahrscheinlichkeitsgesetzz 
@(©)d® verteilt ist, und nimmt man zusätzlich &(©) =const. an, so ergibt sich die Bayesscha 
Wahrscheinlichkeitsverteilung (‚‚inverse probability‘) wz(0|7T)d© für ©, diesichim Gegensatz 
zu II und III auf alle ©, aber nur das eine beobachtete 7 bezieht: ee rn 


ES (T|o SEEN 1. 
wz(9|T) ei TO Be a a ae (8). 
Sseriway 
Eingehendere Darstellung und Begründung dieser verschiedenen Methoden (vgl. [2a, b]) ie 


zeigen, daß die Angabe der Mutungsgrenzen (II) eine voraussetzungs- und einwandfreie, be- 
friedigende Lösung. des Rückschlußproblems liefert und somit vor den anderen den Vorzug 
verdient. Hingegen beantwortet die Bayessche Lösung (IV) eine wesentlich andersartige und 
kompliziertere Fragestellung und ist an die Gültigkeit der Annahme einer konstanten Anfangs- 
verteilung (©) =const. gebunden. Iund IV sind bekanntlich durch die Eigenschaft verknüpft, 

daß das Dichtemittel D, der Bayesschen Verteilung Gl. (8) sich mit der Maximum-likelihood- 
Schätzung Gl. (2) de Jedoch zerreißt diese Bindung, sobald der Parameter ©: einer nicht-- 
linearen Transformation. O=g(£) unterworfen wird (vgl. [2a], S.8). Nur wenn für den trans- 
formierten Gesamtheitsparameter & die entsprechend transformierte Anfangsverteilungsdichte 


5) = wlgk) N DE EN ER IE E 


als konstant angenommen wird, — was jedoch der Annahme &(®) =const. widerspricht, ein 
Umstand, der von R. A. Fisher [1a] (vgl. [2a], S. 7) gegen die Berechtigung der Bayesschen 
Voraussetzung ins Feld geführt wird —, decken sich Dichtemittel der entsprechenden Bayes- 
Verteilung und Maximum-likelihood- Schätzung. 


Es taucht nun die naheliegende Frage auf: „„Besteht auch zwischen den Lösungen I und II 
bzw. — was dasselbe bedeutet — I und III ein ähnlicher Zusammenhang?“ Es wäre zunächst 
zu vermuten, daß irgendeine der wichtigsten Lagekonstanten der Rückschlußverteilung, etwa 
der dichteste Wert oder der Erwartungswert oder der Zentralwert von Gl. (7) mit der Maximum- 
likelihood-Schätzung Gl. (2) übereinstimmt. Im F olgenden wird nun bewiesen, daß der dichteste 
Wert der Mutungswahrscheinlichkeitsverteilung Gl. (7) i. a. nicht mit der Maximum-likelihood- 
Schätzung zusammenfallen muß. 


Unterwerfen wir nämlich den Gesamtheitsparameter © und den Stichprobenparameter T 


°) vgl. [2b], 8.8, 


; Er ‚Partielle Ditferentiation von Gl. ad) nach & 


r 2 
BR a 


een 


LITE 


EN elle alle für die Durchführung der folgenden Reha En, a. Differentiation uneigentlicher 
se Integrale unter dem Integralzeichen — erforderlichen guten Eigenschaften haben mögen, so 
lautet die Dichte der transformierten Wahrscheinlichkeitsverteilung = neuen StichproBedene 
gi parameters 2% bei festem 'Gesamtheitsparameter. repet | 


- 


Me Br. Feld) (0. la) 


Y% 


Die nach Gl. (7) Sachen Rückschlufiyerteilunsdichte von “e aut mi des beobachteten 
gi Wertes 2 Air Eerehen durch: 


fe. . N 


d 
sea - [93 „400 20. | fan | 


5 ACH 9=4(&) 


feld) _dpk) dad) 2f)|O) 
A de t 0° 


9=9(&)- 


und Nullsetzen dieser Ableitung liefert nach GI. (3) für £ die Maximum-likelihood- Schätzung, | 


| eh URla) Fee ae 
die sich auch direkt durch die entsprechende Transformation von Gl. (2) ergibt. Daß Gl. (3) 


tatsächlich ein Maximum liefert, wird durch die wegen Gl. (3,) und (4,) geltenden Gleichungen 
- „ bzw. Ungleichungen ! 


a or füralle 17 Saba ee er ee, 
t= L(z) j 
o2kf(z|d) | _ dp), (dal)? +f(p@)|9) zZ 
SIDE ala > (742) ERFURTER | El) a Se erg (4;) 
t=L(Q t=L() @=L(pı2)) 


gewährleistet. Um dasDichtemittel von w (©|2) zu bestimmen, muß man GI. (7 partiell nach & 
differenzieren: 


2 om (£|z) N (| a) + ae) a 
y 9=g(8) 


Setzen wir nun voraus, daß die ursprünglichen Parameter © und 7 die fragliche Eigenschaft 
besaßen, daß nämlich das Dichtemittel von il (@|T) mit der Schätzung Gl. (2) zusammenfällt, 


“d.h. daß mit einem bestimmten I 


7 
a ET ET Re FA (13)), 
0=L(T) i 
02!m(OIT 
ee Be et, 
®=L(T) 


3) Esist durchaus zulässig, speziell »(2) = 2 anzunehmen, d.h. nur den Gesamtheitsparameter © zu transformieren. 
Jedoch bietet der allgemeinere Fall bei passender Wahl der beiden Transformationen g(£) und »(2) größere Anschaulichkeit. 


dy= BIER: IP @)- un . ee De = 


T 


e En nö 


| Br a 
außer für eventuelle diskrete Nullstellen von en, da ja w(@|7) 


also gewiß positiv ist. Gl. (8) liefert demnach nicht das 
Übereinstimmung des Dichtemittels der 


emittel yon Gl, st le die 
Schätzung ursprünglich vorhanden, so geht sie bei geeigneten Traßkfpumatigmen, der-i m 


Parameter verloren, womit die Behauptung bewiesen ist. 


Ber ee Re 
Ze 5 jf : j a ee = Br ur. vr 
1) -de > 5 ER ER 
. Yan j $ SC aRE < 
mal verteilte Gesamtheit der x-Werte hat den Zentralwert ©. Entnimmt man dieser Gesamtheit Zurück- 
En Stichproben vom Umfange n=3, so ist der Zentralwert 7 der Te Deksenstich ach dran Pie 
scheinlichkeitsgesetz BE n»4h 4 er 


Beispiel: Die nach dem Gesetz 


a ‘ ah 


; T—0 
< 9)" | = 5 
LE 1- | fe *"ay) jar 
20Y2r u : 


verteilt. Die Maximum-likelihood-Schätzung von Oist O=T,da Fr 


T—9 78 
o/(T|®) 3 er 5 \ RE mu -5 
CHI) I EA N Bo 20: ; 
Ra or ni fe 6) I 2 (f u)|++ e e ""dy 
DeTE 0 


für = T verschwindet und 
| &@f(T|®) | 3(r + 4) 


09: BET: 
o=T Y2x 


ist. Die Rückschlußverteilung von © ergibt sich nach G1.(7) zu: 


<0 


Pr 


ma Ka 2.400 22 


7-8 

_-9i ES 

Som“ e elf “a)|=nzıen, 
0 


deckt sich also formal mit der Stichprobenverteilung /(7'|®), und mithin stimmt das Dichtemittel von tw(@]| 7) 
mit der Maximum-likelihood-Schätzung überein. Betrachten wir nun statt © und 7 die dritten Potenzen der 
Medianwerte ®): 


w(O|T) = 


I = h(9)= O8, z=q(T)=KN)=T'%, 
so ist z in Stichproben vom Umfange 3 mit der Wahrscheinlichkeitsdichte 


=T-VE 
WE .Vr-VENR Derset 
J 18) „1: VD) 123 tea .e se 1-5 e ig 
353 20 „8 yon 
verteilt; 
fe ___1 ‚ardzlo) 
oL 9: PX2) 
0=YF 


verschwindet für &=z; das beobachtete z liefert daher die Maximum-likelihood-Schätzung =z. Die trans- 
formierte Rückschlußverteilung von £ ist gegeben durch 


iw(£]z) AAN 
: 3.28 3.68 


*) Die Transformation des Stichprobenparameters 7 ist nicht erforderlich und kann als überflüssi; 
ebensogut unterbleiben, ee 


a 
r 
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ihre partielle Ableitung nach £ durch: 


Amel) _ _2:mWElle , 1 ‚ömtaiyz) 
zu WEB LE 
| e-VE 
Also ist für ö=z : 
= 
= et ee 2 
F 30:28. Yon 
=2z 


das Dichtemittel von i(£]z) liegt demnach nicht in € =z. 


Das Parameterpaar des angeführten Beispiels, Zentralwert einer normal verteilten Gesamt- 
heit und Zentralwert in Stichproben des Umfangs n=3, besitzt die die obigen Rechnungen 
erheblich vereinfachende besondere Eigenschaft, daß die Rückschlußverteilung m(@|7) mit 
der Stichprobenverteilung f(T|®) formal übereinstimmt. Diese wichtige Eigenschaft, die 
natürlich nicht nur das Zusammenfallen des Dichtemittels von w(@]7) mit der Likelihood- 
Schätzung Gl. (2), sondern die Übereinstimmung aller auf dem gleichen Wege aus w(0|T) 
und f(T|@) gewonnenen Größen umschließt, besitzen allgemein alle diejenigen Parameterpaare 
9, T, bei welchen die 7-Verteilung f(T|©) nur vom Absolutbetrag der Differenz |7—0©| ab- 
hängt, also von der Form 


HEN ENT II re 
ist. Aus Gl. (14) folgt nämlich laut Gl. (7) mittels einfacher Integraltransformation 
[are 9) [d f 
m (|T) — ee! =. [FO | 0 9a — | 9az-sr-op), 
‘ DE %=(t-0)% (T-0): 
WROP PS FU Be FATTON a er (15): 


Infolge der Symmetrie der Verteilung Gl. (15) stimmt bei den Parameterpaaren dieser besonderen 
Klasse auch der Erwartungswert M, der Rückschlußverteilung, 


M,=/[9:n(0|Mde=T+Js-fa)de=T......:. (16), 


mit dem beobachteten Stichprobenparameterwert 7 überein. Nur wenn die symmetrische 
Funktion f(x) (welche in <—=0 jedenfalls einExtremum aufweisen muß) in x=0 ein Maximum 
besitzt, d.h. wenn die Maximum-likelihood-Schätzung 


ö=L(N=T 
1° BE 


oy2zr 
der Fall ist), fällt dann demnach auch der Erwartungswert der O-Rückschlußverteilung mit 
dieser Schätzung zusammen, 


lautet (wie es beim obigen Beispiel sowie bei Normalverteilung f(7|®) = 


Der genannten speziellen Klasse gehört allgemeiner noch das vom Gesamtheitszentralwert © einer normal 
verteilten Größe und vom Stichprobenzentralwert 7’ in Stichproben ungeraden Umfangs n gebildete Parameter- 
paar an, da die Verteilung des Zentralwerts 7’ in Stichproben vom ungeraden Umfang n, 

T—09 
r FE=0 y® 27 n—1 
! 20: 20: 2 
ee an ir te sa er .dy i 


2) 27-0 12% 7% 


nur von (T—9)? abhängt. Alle nur aus dieser Eigenschaft für n—=3 gefolgerten Eigenschaften der zugehörigen 
Rückschlußverteilung gelten daher auch allgemein für beliebiges ungerades n. 

Im allgemeinen deckt sich der Erwartungswert M, der Rückschlußverteilung mw (©|7) 
nicht mit der Maximum-likelihood-Schätzung Gl. (2), wie man ebenfalls am obigen Beispiel 
erkennt: 


m Geppert, cin eo hung vn cr 


Der Erwartungswert Mg der im obigen ae ereananer Bene Rückschlußyerteilung i(£]e) £ 
ergibt sich airktela Seenforneken der Integrationsvariablen und hc en Ze ee: zu: u 


12 -  _e-m a = Se Bi 
ER ne En | 12. )| 3e 7 € 
2oy2r 2% 7 


| | | N ' Rt 

2 er, ee 

3 3743712419: e "in l elf en er) Re =2 

SPILE fi ° l no® : $ FRE Ei 
© x i k 

a Yu’ . 
ui Tr, 2 20: aan BEA A 
fs x + Ta) e 20 . Inst . (J- ) | dx = >, (3 0® + T?a 03 e7 5 


0 
M:=T-(M+30)-2+30-Ye+z, 
fällt also für 2-0 nicht mit der Schätzung {=z zusammen. 

Die Frage, ob. der Zentralwert Z, der Rückschlußverteilung mit O=L(T) zusammen- 


fallen muß, läßt sich nicht mit dem oben verwendeten Hilfsmittel der Parametertransformation 
beantworten; denn übt man eine Transformation Gl. (10) aus, so ergibt sich der Zentralwert Z; 


der transformierten Verteilung Gl. (7) wegen 


. h(Zo) h (Ze) B Ze 
[Beia«= [vWolre)- Dur- [» (0|T)40=; 


offenbar, indem man den ursprünglichen Zentralwert Z, in derselben Weise transformiert: 
Z;=h(Ze). 

Deckt sich also der ureptinpliche Zentralwert Ze mit der Schätzung L(T), so bleibt diese Eigen- 

schaft auch bei Transformation der Parameter erhalten. 


In dem oben benutzten Beispiel stimmt der Zentralwert Zg von w(®]| T)tatsächlich mit der Schätzung 
nach dem Maximum-likelihood-Prinzip überein, denn es ist: 


m Yu ö 9—T 
te=T) 2 ir 
SL Pr De RT I er \o- 
20 Y2n tal 
00, 00 0 
SQ z 
3 x: y? 2 
ml 2 feel ee 
20 Yan “ n0® 2oY2ar 2 6 2° 
u 0 


Folgendes einfache Beispiel zeigt jedoch direkt, daß die Maximum-likelihood-Schätzung 
Gl. (2) nicht notwendig den Zentralwert der Rückschlußverteilung GI. (7) liefert. 


Es sei O=0? die Streuung einer nach dem Gesetz 
(z—a)? 
1 FB, 
.e 
oy2r 

um den Mittelwert « normal verteilten &-Gesamtheit. Die in der üblichen Weise in bezug auf den Stichproben- 

mittelwert berechnete Streuung T' in Stichproben vom Umfange n =3, die dieser Gesamtheit entnommen werden 
folgen nach Helmert der Verteilung 


dx 


3.0.28 
KTIOdT=-,5.e dt, 8>0,T20 
Wegen 
er 
of(TI® i 29 
Li Er (37 —20) 
ist die „‚beste‘‘ Schätzung von © auf Grund von T: 
rm 
2 
Er he 
en ee 
I * 


er 


x 


Pu: 
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EUER RER EEE 
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Die Rückschlußverteilung ergibt sich nach c1. (7) als : a Be Er Ze 
Wr 
eh Mr ar f | ne 2 38 i t } n ; 
er RER N Be ee Re ze 
.das-von 0 bis zur „‚besten‘‘ Schätzung L(T') erstreckte Integral derselben 
| ner ST. ma Er 


SEE Sg > ER 37 je Ko 
Y; x / x . Zar a 2 
2 R 3 a a 3 Bee —1 |] 
Sramo=2.2. | ro-tr. Er 


Re OF EE © RE 
A R f 12% Er s ’ j ; ß 
ist. verschieden von >> u.h, O= 5 T ist nicht der. Zentralwert der Rückschlußverteilung. 
Übrigens liefert dieser Fall auch ein besonders einfaches Beispiel für das mögliche Abweichen des Dichte- 
mittels der Rückschlußverteilung, welches hier 7 = 2: T lautet, von der „besten“ Schätzung © -> T. 


Trotzdem wurde dem obigen umständlicheren Beweisgang der Vorzug gegeben, da er gleichzeitig Einblick in das 
Verhaltender Lagekonstanten der Rückschlußverteilung gegenüber Parametertransformationengewährt. Übrigens 


läßt sich der Erwartungswert der Rückschlußverteilung in diesem Spezialfall nicht ohne weiteres bestimmen, 
oo 


da das uneigentliche Integral [ .o-w (0|T)d® offenbar am unteren Ende nicht konvergiert. 


5 


Es ist somit gezeigt, daß weder das Dichtemittelnoch der Erwartungswert noch 
der Zentralwert der Rückschlußverteilung mit der „besten“ Schätzung nach dem 


Maximum-likelihood-Prinzip zusammenfallen müssen. 
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KLEINE MITTEILUNGEN : 


Über das Verhalten eines angefachten schwin- 
gungsfähigen Systems mit einem Freiheits- 
grad, dessen Dämpfung dem Quadrat der Ge- 
schwindigkeit proportional ist. 

l. Fragestellung und Methode. Die 
Bewegungsgleichung des Systems laute 
+nasgne+ixe=xsint (n>0, A>0, x>0), 
wobei YA das Verhältnis von Eigenfrequenz zu Kraft- 
frequenz bedeutet und nach Lösungen mit der Periode 
2n gefragt ist. 

Es werde einschränkend vorausgesetzt, daß 

1. x so klein sei, daß höhere Potenzen gegenüber der 

ersten hinreichend zurücktreten, 

2. &(t) für 0 St < 27 genau zwei Nullstellen habe. 
Es dürfen dann. diese beiden Nullstellen nach it = 0 
und t= zgelegt werden, so daß die Bewegungsgleichung 


„unter Einführung eines nachträglich zu bestimmenden 


Parameters o in der Form geschrieben werden kann 
z+n&:sgnsint+Ax=xsin(t+o). .(l). 
Eine Grundlage für die Behandlung von Gl. {1) bietet 
die Integralgleichungstheorie von E. Schmidt!). 
Die zu Gl.(1) gehörige Integralgleichung (Integral- 
potenzreihe) kann mit &= u ersichtlich in der Form 
27 
ul) + [IK (9) u (0) + K,(t, 5) (u (s))?]ds 
= +x(cos(t +0) —coso)=0 
ı) E. Schmidt, Math. Ann. 65, 370—399, 1908. 


geschrieben werden und es sind nun Nachbarlösungen 
zu der zu x = 0 gehörenden Lösung  (t) = 0 gesucht. 


Diese Nachbarlösungen werden nachE. Schmidt 


außerhalb der Eigenwerte der«zugehörigen linearen 


Gleichung ä+Au=0 bzw. @&+Ax=0 in erster 
Näherung durch Wegstreichen der nichtlinearen Glieder 
gewonnen, und jede folgende Näherung xm+1 ergibt 
sich aus derjenigen inhomogenen linearen Gleichung, 
die durch Einsetzen der vorhergehenden Näherung xm 
in die nichtlinearen Glieder entsteht. Es kommt dies 
auf eine Potenzreihenentwicklung nach 7 hinaus. 

Ist hingegen A ein Eigenwert, so werden in der Fou- 
rierentwicklung der rechten Seite von (1) die zu dem 
betreffenden Eigenwert gehörenden Fourierglieder in 
der inhomogenen linearisierten Gleichung unterdrückt 
und nachträglich mit unbestimmten Koeffizienten zur 
Ausgangslösung hinzugefügt. ä 

Diese hat demnach für A = 1 die Form 


Acost+ Bsint 
und für ment (M=2,) 


sin (t+ 0) + Acosnt + Bsin nt. 


Die unbestimmten Koeffizienten A und B werden 
durch die Orthogonalitätsbedingungen- (Verzweigungs- 
gleichungen) festgelegt, die die Auflösbarkeit der Glei- 
chung für x, 


&, +n?a, = xsin (+0) —nalsgnsint 


N = - NER 
a 


NS 
“. 


zum Ausdruck bringen, also 


2n x 
a #?sgn sin t)-cosnt dt=0 


2n 
f-..)-sinnt dt=0 
6 % 
Der Phasenwinkelo ergibt sich schließlich aus 
&,(0)=0. 
"Die höheren Nähe; en für & (t) werden durch sinn- 
gemäße Wiederholung des Verfahrens gewonnen, 
Es interessieren hier nur die ungeraden Eigenwerte, 
da für A=4,16,.. . die Gl. (2) auf A= B = 0 führen, 
2.Ergebnisse, Ergänzt man, um auch die Um- 
bung des Eigenwerts zu erfassen, Gl. (2) noch durch 
heits- und Rückstellkraft, so erhält man für} 1 


mit, = Acost+ Bsint 
3n(A4-—1)\2 


em en 
Ben” 
a ea) 


67% 
und hieraus für das Resonanzmaximum bei A=1 
unter Hinzufügung der aus der zweiten Näherung hin- 
zutretenden Oberschwingungen 


Ehen EL En A 7 
= —, grnrcost— Inn% PETTeNGEr 


y=B3,5... 


(3) 


+:4’ cost + B’sin )| + ol«®), 


wobei A’ und B’ von n und x nicht mehr abhängen. 

Das negative Zeichen von A in (3) bringt zum 
Ausdruck, daß gemäß (1) x (t) für 0 <t <x positiv 
sein muß. 

Für die „Verstimmung‘‘ A o, um die die Frequenz 
@,. der anfachenden Kraft von der Eigenfrequenz @; 
des Systems abweichen muß, damit das Amplituden- 
quadrat auf die Hälfte des Resonanzwertes sinkt, folgt 
aus (3) 

Aw n%* 
FA ke Re 


3 
Für }=9 wirdo = 9y7 und 
% 
2 = — 7 (cos 1 + 0,042 cos 31) + 0 () 
und für = 250 = 1e und 


% 
© = — 54 (cos t + 0,0038 cos 51) + O (x). 


Wesentlich für die quadratisch gedämpfte Schwingung 
ist also 

I. das Auftreten von Oberwellen ungerader Ord- 
nung, wobei die n-te Harmonische besonders her- 
vortritt, wenn das Verhältnis von Erregerfrequenz 
zu Eigenfrequenz 1/n ist, 

2. das Anwachsen der Amplitude und das Absinken 
der Resonanzschärfe mit der Quadratwurzel der 
Erregerkraft. 

Die Beobachtung dieser Erscheinungen dürfte infolge 
des Umstandes, daß Strömungswiderstände bei kleinen 
Geschwindigkeiten proportional der ersten Geschwin- 
a anwachsen, wohl nur qualitativ möglich 
sein ?), 

Die Konvergenzgrenze des benützten Verfahrens 
könnte an Hand einer eingehenderen Betrachtung der 
angegebenen Integralgleichung abgeschätzt werden. 
Insbesondere gibt es eine Zahl cı > 0,so daß die Reihe 
mit dem Anfangsglied (3) für n*<c, und etwa 
0.<A<3 konvergiert,. 


°?) Handb. d. Exp. Phys. Bd.17 Aufsatz Martin. 


A 


3. Verwandte Formen des hr es 


A Ist die dissipative Kraft u; 
dritten Potenz der Gecshwindigkeit, proportional, gilt r- 


also ve 2: 

i+n®+ıiz=xsint, : z 
so wächst für A=1 die Resonanzamplitude mit 
Kubikwurzel der Erregerkraft an: 


/E 3 


Weitere Einzelheiten bieten hier kein Interesse. 

Ist die dissipative Kraft i 
(Reibung fester Körper), so setzt bekanntlich erst ober- 
halb eines gewissen Schwellenwertes der Erregerkraft 
die Schwingung ein, erreicht dann aber rechnerisch 
sofort eine unendliche Resonanzspitze. Auch hier treten 
Oberwellen auf. 


Darmstadt. G. Plato. 
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Zurückführung der Wärmespannungen in 
einem elastischen Körper auf ein Knick- 
Biegungsproblem. 
’ Einleitung. 
Wir bezeichnen mit «, v, w die Komponenten des 
Verschiebungsvektors v, mit e=divv die kubische 
Ausdehnung, mit = #(x, y,z) das gegebene statio- 
näre Temperaturfeld, dem der elastische Körper aus- 
gesetztseinsoll. Dann gelten bekanntlich diefolgenden 
Differentialgleichungen 
ar Dr PR 
m—2 0: m-—2 dr (x, m Kon- 

An mr 00 m niet. di 
m—2 d0y: m—2 dy Körpers) 

A 2(m+1)a 08 m de 

= ma . m—-2% fe 
und ferner an der Oberfläche die Randbedingungen, 
die das Verschwinden der Normalspannung und der 
Schubspannung aussagen: 


ou e m-+1 
24 E - N x) 
ou 
+4 (5 + n cos (n, y) 


B ow 
+4 Pa) Pr cos (N, 2) = 0 


E£ ee 
@ Ep + 3 y, cos (N, &) 


ov e m+l 
+2@ [e Tr a ya. 0)cosn, Y) (2) 


o9v dw 
+@ 82 39, cos (n,2) = 0 


o9w du 
@ ( > I an cos (N, &) 
0w Ov 


+4@G E- 4 & cos (n ‚y) 


fu e ml 
+2@ pr +33 cos (n,2)—=0 


(r innere Normale, G konstant). 


Fu 
das Glied 37 Ep durch m 


ae = 
du 
ee Bee te 2 


re a = ou = 


= re =7,ugedu 


n= (cos (n, x), cos (n, Y), cos (n, z)) 


oo. 


du 
lm +) cos (m, y) 


( 
Een. 
+(&- 


)' 
Tr = + dy (cos (n,9) 


m 


— > + 08 (n, y) rotz » — cos (n, 2) rotyd. 
=2n- - grad u +[n x rotole. 


=) cos (n, 2) 


die Randbedingungen in der Form 


z|n- -gradv + Sarnen (2a) 


+[n x rot] = 435 R 
geschrieben werden können. 

Es wird im folgenden eine vereinfachte Schreibweise 
der Grundgleichungen angegeben, durch welche die 
Zweiteilung des Verschiebungsvektors nahegelegt wird, 
auf der der auszusprechende Satz beruht. 


Beweis des Satzes. 
Ersetzt man in der ersten Grundgleichung gemäß der 


Qu 20 | du 
= trayr% 
Pe kn 
& un EB erhält 


' Bedeutung 


man: 
de 2m—-2. Ru 0% Ru Ow 


O5 ma tape duoy ' Oz due 
2(m+1)a0% 
- m—2 9z 
oder : - 
2 m — 2 ae ee) m+1 08 
2m 1)\ 94 02 m—1*9x° 


ER kann man die Grundgleichungen in der Form 
schreiben 


—2 m+]1 
grade — She (rott)= 1% grad® (la), 


ee 
oder, wenn man 


en rtv =q=D .'. (3) 


setzt, auch so 
ee pren-0 ee th} 


Es liegt nun nahe, eine > et (3) entsprechende 
Zerl des Verschiebungsvektors selbst zu suchen. 
Dazu braucht man nur zu setzen: 


v=n, +% 


rl | oona;,0 | 


wobei das Integral über den gesamten Bereich des Kör- 
pers erstreckt ist und P den Integrationspunkt, Q den 


bedsatehs findet man unter Benutzung der 
 Poissonschen Formel: 


ne A 


en 4 a 1 
m 


also wegen (3) und (4): FE 
div® = divv —divv, —n, 
Da das Vektorfeld p, ein Eoteatıag besitzt, ı so ist 


rotp, = (0, also. 


rot® = rot —rotd, = rt =g=QD. 


Die Gleichungen (1b) können daher als die Grund- 


gleichungen im Falle eines’ konstanten Temperatur- 
feldes gedeutet werden. Es fragt sich, wie nun die zu- 


_ gehörigen Randb en aussehen. Einsetzung 
yon Gl. (4) in (2) liefert: 
divv,+H 1 
2 hr (grad vd, + grad®) + ne on] 


+fhxD]= 
oder nach kurzer RE 


2 n- maß +, 2 | + x rot] 


- 201 sdn—ı-gradv,]. 


Man sieht demnach: Das Vektorfeld ® stellt die 
Verschiebung für den Fall konstanter Temperatur Null 
dar, wenn auf die Oberfläche des Körpers die Ober- 
flächenkraft 


260 1m jadn—n- grad d; | 


wirkt. Diesist aber ein reines Knick-Biegungsproblem, 
das nach den in der Mechanik üblichen Methoden be- 
handelt werden kann. Unser Satz lautet demnach: 
Die Verschiebungen in einem Kör- 
per, der von der Temperatur Null 
auf die Temperatur %(z,y,2) erwärmt 
wordenist, lassen sich aus zwei An- 
teilen zusammensetzen. Der erste 
Anteilergibtsichnach dem Gesetz: 


m+1 & 
-m—1 [Here 


Der zweite Anteil ® stellt die bei der 
Temperatur Null eintretende Ver- 
schiebung dar, welche durch die 
Oberflächenkraft 


2G 


. vd, 


m+]1 
ms dn—n-gado, 


erzeugt wird, und ergibt ein reines. 


Knick-Biegungsproblem, welches 
nach den in der Mechanik üblichen 
Methoden behandelt werden kann. 

Zur Erläuterung betrachten wir als Beispiel 
den Falldes Hohlzylinders. Der Vektorv, 
wird proportional dem Kraftvektor eines Schwere- 
feldes, dessen Massendichte dem Temperaturfeld 9 pro- 
portional im Hohlzylinder verteilt ist. Dieser Anteil 
kann durch einfache Quadratur in bekannter Weise 
ermittelt werden. An der Oberfläche hat man weiter- 
hin senkrecht zu Deck- und Mantelflächen den Vektor 

2G an j%®n anzubringen und davon den Vektor 
2 en n* Er p, abzuziehen. Dabei ist bei letzteren je- 
weils der innere Randwert zu wählen, da ja der stetige 
Anschluß nach dem Innern zu gewahrt bleiben muß, 


Berlin-Frohnau. W.Schmeidler. 


rg die Entfernung zwischen en Di 


ET 
+ 


u  } an | [Jeanne 2 ä 


Br 


dauer und Rückstellkraft. ET 
Ist bei einer Schwingung von der Form —- = 


_ Ya)—V(a) die Rückstellkraft nichtlinear, das Po- 
'tentia] also nicht von der Form 1/2 cx?, so kann die 
"Schwingungsdauer T in Abhängigkeit von der Ampli- 


tude a inkeinem endlichen a-Intervall konstant sein, 
Denn für V, (x) =. Ja}? ist 


e) 12 
RE 5 
= 


und für V, (x) = eat +42+ + mean) ist. 


RR Ze il, 39 
T,(a) 2. Yrlı— za a (Uu-8e) alt: \e) 


mit endlichem bzw. bei positivem c, sogar unendlich 
großem Konvergenzradius. ü 

Soll T(a) für O<a<a, konstant sein, so muß 
daher das als ganz-rational (oder auch analytisch) vor- 
ausgesetzte Potential notwendig die Form 1/2 c@* 
haben. | 

Weiter zeigt die durch Variablentransformation 
leicht zu gewinnende Darstellung 


1 
7 =YEm:2a | yay 
0 


kr 
T,(@)=2 "an (2) 


YV (a)— V(a—ap?) 


die Differenzierbarkeit und damit das reguläre Ver- 
halten und die Fortsetzbarkeit von T7(a) auf der 
positiv-reellen a-Achse. Hieraus folgt, daß bei nicht- 
quadratischem Potential auch in keinem Intervall 
der Form a, <a <a, die Funktion T(a) konstant sein 
kann. E 

Gl.(2) kann u. U. dazu dienen, zu gegebenem T7'(«) 
das zugehörige V(x) zu finden. 


Darmstadt/Hessen. G. Plato. 


Die Hertzsche Gleichung für elliptische 
Druckflächen. 

Daß die Hertzschen Gleichungen für den zwei- 
dimensionalen und rotationssymmetrischen Fallrichtig 
sind, ist leicht nachzuweisen. Anders verhält es sich 
für den allgemeinen Falleiner elliptischen Druckfläche. 
Hierfür wird der Beweis nach Hertz mit Hilfe des 
Potentials eines Ellipsoides durchgeführt. In vorlie- 
gender Mitteilung möchte ich einen elementaren Beweis 
bringen, den ich in meinen Vorlesungen über Elastizi- 
tätslehre in Dresden mehrfach vortrug. Sollte er be- 
kannt sein, so bitte ich, mir die Literaturstelle mitzu- 
teilen. (Im übrigen hoffe ich, daß die leichte geistige 
Kost den eiligen Lesern munden wird.) 

Die Oberfläche des Halbraumes z > 0 sei innerhalb 
einer Ellipse mit den Halbachsen «a und b belastet durch 
den Druck 


1/ za /(Yyı2 
Pr li) -(E) - .- ©. 
Die Druckkraft wird hierbei 


2 
P= q RT 1 Ko (2). 


Zu beweisen ist, daß die Verschiebung der Oberflä- 
chenpunkte in z-Richtung innerhalb der Ellipse 
vn tn tun... (3) 
wird. Außerdem sind die Koeffizietenn a,, a, und a, 
zu bestimmen. Der Beweis wird ähnlich wie beim ro- 
tationssymmetrischen Fall durchgeführt. (Siehe z. B. 
n I ppl, Vorlesungen über technische Mechanik, 
AV): 
Ist der Halbraum durch eine Einzelkraft P belastet, 


ber den Zusammenhang von Schwingungs- 8 


. Im Bild 1 ist ein solcher Schnitt eingezei 


Für eine beliebige Belastung p (z, y) wir 


2n © 
1 1— 2 (P 
u. [rd | Zrapar 
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Enz 
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+ er. 

Das Integral [pdr kann man sich leicht vorste- 

Br a 

len als Schnittfläche des Belastungshügels p, wobei der BR. 
Schnitt durch den Punkt A unter dem Winkelggeht. 


2 Du a a 


Bild 1. Beliebiger Belastungshügel mit Schnittfläche 
unter dem Winkel, 


Nunmehr gehen wir zu dem Belastungshügel nach 
+o 


Gl. (1) über. Um das Integral f pdr zu berechnen, 


vergleichen wir die belastete elliptische Fläche mit einer ““ 
Kreisfläche vom Halbmesser a. In den Bildern 2a und 
2b sind beide Flächen dargestellt. Die Kreisfläche ent- 
steht aus der Ellipsenfläche durch Streckung in y-Rich- 
tung im Verhältnis a:b. Jedem Punkt der Ellipsen- 
fläche entspricht ein Punkt der Kreisfläche; die Be- 
lastung in den entsprechenden Punkten sei gleich groß. 
Alle Elemente der Kreisfläche bezeichnen wirin gleicher 
Weise wie die Elemente der Ellipsenfläche, nur mit ge- 
strichenen Buchstaben. 
Hierbei ist , 


und 


3 x\2 y'\a 
pP = Bm ]/ 1 (7) (5) Sl, 
Nunmehr wählen wir innerhalb der Ellipse einen 
Punkt A, dessen Verschiebung w wir bestimmen wollen. 


Hierzu führen wir im Grundriß eine Schnittgerade 
durch A unter dem Winkelp und berechnen erst das 
+ © 


in Gl. (5) vorkommende Integral [ pdr., 


stläche der Punkt A x 5 
Die entsp Er, Sehntigend ern unter 


mel p. 


ierbei ist 
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Bild 2a und 2b. Belastete Ellipsenfläche und Kreisfläche. 


Der Abstand vom er pn bis zur Geraden 
durch 4A’ wird 
x sing’ — y’ cos’. 


f Die halbe Sehne hat die Länge } 


er TEE \2 

Pe: r [7 1— (7 sinp —, 0sP) . 

£ -— Die Höhe in der Mitte der Schnittfläche wird 

a x’ ; , Yy 12 

Dsr 1— (7 sinp cos) - 

Die Schnittfläche wird also 

E14 En y' Ns 

5 Pinax ® Er sing, — „cosp) | 


Die entsprechende Schnittfläche des Belastungs- 
“ hügels über der elliptischen Fläche erhält man hieraus 


durch Multiplikation mit er Also 
+ ©. 


dr=> 1 RR: — 4 cos Be 
DR Pax @ er ing p c08 9 
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2 =7 5 Yundführen an Stelle von p ‚die Neradalrueg 
ein. as G1. (8) folgt 

Be BAT 
cos? p PT a cp rl 


Weiter wird cosp nach Gl. 2 ersetzt. Mit diesen 
Umformungen_ wird 


De ee * 2 x Ba Be 
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Das ge zerlegen wir in vier Teilintegrale und 


erhalten 
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y1— k2sin?p’ 


K und Esind die vollständigen elliptischen ge 
1. und 2. Gattung. 


Wir erhalten das Endergebnis 


1— x \ 
vg "Pad |E- al | e 
(18). 
 al@lK+ 2) 
(20b) Schlewecke ü. Bockenem. C.Weber 


Um das Integral auszurechnen, Reigen wir « = 2, 4 
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L.J. Comrie, M. A., Ph.D., Chambers’s four- 
figure mathematical tables. Edinburg 
a London 1947 bei W. u. R.Chambers Ltd. Preis: 
in blauem Umschlag 5, in Leinwand geb. 6 Sh. 
Ein Versuch, eine vierstellige Tafel zu schaffen, die 
- sowohl für. die Oberklassen der höheren Schulen, wie 
für die Hochschulen und die industrielle Praxis geeignet 
ist. Die Tafel gibt Logarithmen, Antilogarithmen, 
Kreis-, Hyperbel- ünd nentialfunktionen; Qua- 
drate, Quadratwurzeln und Kehrwerte aller dreistelli- 
gen Zahlen (bei allen Kommastellungen); Ga u ßsche 
Verteilungsfunktion und Wahrscheinlichkeits-Integral; 
dazu viele Hinweise und Erläuterungen, 
Bemerkenswert an dem Buch sind die Einrichtung 
und die Ausstattung, also die Tafeltechnik, die auf die 
roße Erfahrung Comries in der Herausgabe von 
Zenlesarlne ründetist. Das Buch enthält 64 Seiten, 
es ist 195 mm breit und 255 mm hoch. Die Form ist 
also viel größer als sonst bei vierstelligen Tafeln. Es 
ist nicht in einen steifen, sondern in einen biegsamen 
blauen Umschlag gebunden. Auf den Innenseiten des 
Umschlags stehen Proportionalteile, auf seiner Rück- 
seite steht das Inhaltsverzeichnis. Innerhalb einer Tafel 
finden sich weder horizontale, noch vertikale Linien, 
Nur deutlich erkennbare Zwischenräume trennen be- 
nachbarte Zahlen voneinander. Wo die Zahlen h>ri- 
zontal aufeinander folgen, trägt die erste En die 
Überschrift 0, die letzte Spalte aber nicht die Über- 
schrift 9, sondern 10, sie ist also eine Wiederholung der 
ersten Spalte, damit man bequemer die Differenz bilden 
kann. Langsam veränderliche Funktionen werden in 
sog. kritischen Tafeln gegeben, d.h. es wird angegeben, 
zwischen welchen Grenzen der unabhängigen Veränder- 
lichen ein bestimmter Funktionswert gilt. Die Winkel 
werden in dreierlei Einheiten angegeben: 1.in Graden 
und Minuten, 2. in Graden.und Dezimalteilen des Gra- 
des, 3. in Radianten (im Bogenmaß). Bei manchen 
Tafeln werden die Interpolationszuschläge in rechts 
daneben verstreuten Zeilen angegeben, was in engli- 
schen Schulen gebräuchlich zu sein scheint. Wo die 
Zahlen vertikal aufeinander folgen, sind die Differenzen 
in kleinerer Schrift rechts daneben in die Zwischenzeile 
gedruckt. Wenn in einer Zeile vor die Ziffern die ab- 
gespalteten Anfangs-Ziffern der nächsten Zeile zu setzen 
sind, so wird das nicht durch einen Stern angezeigt, 
sondern der größeren Deutlichkeit und Auffälligkeit 
wegen durch schrägen Druck der Ziffern. Der Raum 
ist aufs Außerste ausgenutzt. Jede Stelle, die sonst 
frei bleiben würde, enthält eine nützliche Angabe. 
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Buchhenrecungen — Hngegangene acer En 
BUCHBESPRECHUNGEN 
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Es ist nicht, angängig, hier alle Eigentümlich 


RE L.v.Schrutka hat sich ebenso eı 
schieden für die gleich hohen . 
veränder. 


un licher 
in den Tafeln vonSchrön 
sind. 


ctg 0,5° = tg 
man durch tg 0,5° = 
dieser letzte Wert freilich nicht als 0 


werden, also zwar mit vier Dezimalen, aber 
zwei geltenden Stellen. Es macht keine Schwierigkeit, 


Logarithmen benutzt (3,9409 abziehen statt zu- 
zuzählen). Ohne eines dieser Hilfsmittel wi 
jemand trigono isch rechnen. Aber in diesem sel- 
tenen Falle kann man den großen Multi aus 
einer Tafel der Kehrwerte ablesen, um die Division mit 
einer vielstelligen Zahl zu vermeiden. Im vorli 
Buch würden so die Seiten 20 21 26 bis 29 32 bis 35, 
54, 55 wegfallen, im ganzen 10 Seiten. Comrieselbst | 
weist darauf hin, daß auch die Arcus-Funktionen be- j 
quemer zu finden sind, Se nn von kleinen und 2 
nicht von großen Werten des Tangens oder Cotangens - 3 
eht. Rede steht auch wohl nur ein überliefer- 

treben, an einer starren Einförmigkeit des Tafel- 
bildes festzuhalten. 

Wie man auch zu den Einzelheiten stehen ‚in £ 
jedem Fall wird das vorliegende Buch als Vorbild ein E 
wertvoller Berater für die neuer Tafeln sein. 

Vor allem ist aber das Buch ein nützliches Rechen- 
mittel, an das man sich schnell gewöhnen wird und das 
man dann nicht mehr wird ra wollen. 

Fritz Emde, 
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Bei der Schriftleitung sind folgende Bücher ein- 
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Dr. ing. Ludwig Föppl, (o. Prof. a. d. Techn. Hoch- 
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die rollende Reibung. 448. m. ı3 Bildern. 
München 1947, Leibniz-Verlag bisher R. Oldenbourg- 
Verlag. Preis brosch. 6,— M. 

Dr. Felix Strecker, Die elektrischeSelbst 
erregung mit einer Theorie der ak- 
tiven Netzwerke. 142 8.m.62 Abb. Stuttgart 
1947. Verlag $. Hirzel. Preis brosch. 7,80 M. 

Dr. Fr. Hund, (o. Prof. an der Universität Jena)» 
Einführung iin die theoretische Phy- 
sik. Band II: Theorie der Elektrizität 
und des Magnetismus. Meyers kleine Hand- 
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bücher Bd. 47/48, 297 8. m. 75 Abb. Leipzig 1947. 
a Bibliographisches Institut. Preis Go 
5,80 M. 

D.. Fr. Hund, (o. Prof. an der Universität Jena), 
Einführungin die theoretische Phy-_ 
sik. Band III: Optik. Meyers kleine Handbücher 
Bd. 49/50, 208 $. m. 100 Abb. Leipzig 1947, Verlag: 
Bibliographisches Institut. Preis brosch. 5,80 M. 


Druckfehlerberichtigung: 
In Band 25/27, Heft 7, Seite 189 Gleichung (40) 
muß es heißen: {K}statt K. 
Seite 191 Gleichung (3a) ist ®= w,x statt P=w,'* 
zu Setzen. 
Heft 8/9, Seite 216 in Gleichung (1.2)solles o(1+%). 
statt o(1+%, heißen, 
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